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I. INTRODUCTION 



This dissertation addresses what appears to be a simple problem: “Given a real, 
square matrix A, when does there exist a nonnegative matrix B such that B _1 AB has its 
off-diagonal elements nonnegative?” Based on different bodies of literature, and on the 
application, we could similarly ask when a square linear operator is quasimonotone 
nondecreasing with respect to a cone in the nonnegative orthant, or when a real matrix is 
essentially nonnegative under a nonnegative change of basis. 

This question has its roots in the Perron-Frobenius theory of the early 1900s, and 
as we develop our solution, we trace its history, as well as the history and current state of 
some of its applications. We also discuss similar and related problems to demonstrate how 
this problem fits into the scheme of applied mathematics. 

We present a brief background on cones in Euclidean space and on the theory of 
nonnegative matrices in Chapter n. We use the following convention to present results. 
Theorems (and proofs) from the original source have the author and date after the 
theorem number. If the proof comes from another source, we indicate that prior to 
presenting the theorem. For well-known results, we either indicate the source or provide 
our own proof We indicate the completion of a theorem and proof (or just a theorem for 
results presented without proof) with the □ symbol. 

. Chapter HI motivates this question as an unsolved problem from the theory of 
cone-valued vector Lyapunov functions, which can determine stability in dynamical 
systems. We also present some of our results concerning this application. This question 
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was first posed for the application of vector Lyapunov functions in 1974, and partial 
answers were given for this setting in 1980 and 1995. 

Chapters IV and V present our original solutions to this problem. The problem 
has significantly different solutions depending on whether the spectrum of A is real or 
general, and Chapter IV presents our solution for the real spectrum In Section IV. A we 
show the necessary condition that the first eigenvector xi of A (associated with the 
greatest eigenvalue ?w) be nonnegative, and the sufficient condition that Xi be positive. 

The case where xi is nonnegative but not strictly positive has different solutions 
based on the reducibility of the matrix A. Section IV. B addresses the reduced case, where 
we show sufficient conditions when \\ is a simple eigenvalue of A, and we have a 
sufficient (but not necessary) condition when X\ is not simple. 

In the case where xi is nonnegative and A is irreducible, Section IV. C shows 
that we can use xi to deflate the matrix A and reduce the problem to one of dimension 
n - 1. Here, we also give necessary and sufficient conditions for a cone to exist. 

Chapter V presents our solution for the complex spectrum We show that in the 
presence of a positive first eigenvector, the problem can be reduced to the nonnegative 
inverse eigenvalue problem, a classic unsolved problem from theoretical linear algebra. 

We present the solution only in terms of irreducible quasimonotone matrices, and we show 
why it is sufficient to consider such matrices. 

Chapter VI presents other direct and related applications from control theory and 
dynamics, as well as from other fields of applied mathematics, and we mention similar 
problems which have been solved or discussed previously from these and other fields. 
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Chapter VII summarizes our results, to include the unsolved problems mentioned 
above, and we suggest further research. 

One attractive aspect of our solutions is they are constructive, and hence fairly 
simple, although there is some numerical instability inherent in some constructions. 
Particularly for applications where the existence of a matrix B is sufficient information for 
the analyst, we present collections of useful necessary and sufficient conditions. We begin 
with some background. 



n. CONES AND NONNEGATIVE MATRICES 



This chapter provides background for the main results. Specifically, we discuss 
cones in Euclidean space and present results from the theory of nonnegative matrices. We 
begin with a discussion of cones. 



A. CONES IN EUCLIDEAN SPACE. 

The employment of cones as subsets of R" began in the 1930s, where along with 
convex polyhedra, they were soon found to be useful in the field of linear optimization and 
in the study of linear inequalities (see, for example, Weyl, 1935; and Gol dman and Tucker, 
1956). Dual cones were first employed in 1941 by J. Dieudonne (see Lay, 1982), 
although dual sets to convex sets had been used much earlier (see Helly, 1923). 

The analysis of cones, convex polyhedra, and convex sets in general has become a 
major field of study in its own right, which has its roots in Caratheodory (1907), and 
which can be found, for example, in Sandgren (1954), Karlin ( 1968), Rockafeller (1970), 
Stoer and Witzgall (1970), Berman (1972), and Lay (1982). We need only the most basic 
definitions and results from those fields, and we begin with some definitions. 

If R" is Euclidean n-space with norm || j| and inner product (•,•), a subset 
ScR" is convex if and only if x, y e S implies ax + ( 1 - a)y e S for all a € [0,1]. A 
set K is a cone if and only if XK cK for L > 0. A convex cone satisfies both 
definitions, and we require that K = K , where K is the closure of K, so that our cones 
are closed. (It is a result of Farkas (1901) that a set is a closed, convex cone if and only if 
K = (K*)*, where K* is the dual of K, to be defined shortly.) A cone is pointed if and 
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only if Kn (-K) = {0}, and solid if and only if K° is nonempty, where K° is the interior 
of K An equivalent notion (in Euclidean space) to solid is reproducing, where K is 
reproducing if and only if K - K = R". (It is a result of Krein and Rutman ( 1948) that a 
cone is pointed if and only if its dual is solid.) 

A convex cone K is polyhedral if and only if it is generated by finitely many 
vectors, and if the number of independent vectors is equal to n, then K is simplicial. A 
closed, pointed convex cone is called proper. The result that a proper cone is generated 
by its extremals is a special case of the Krein-Milman Theorem. The only cones we 
consider for our applications are proper, simplicial cones. 

Since a proper, simplicial cone is generated by n independent extremal vectors bj, 
we consider the nonsingular matrix B = [bi, . . ., b n ] e R nxn , and we denote the cone 

D 

generated by the columns of B as K(B). Clearly, K(B) = {x e R n |x = cojbi, Qj> 0}. 

i=l 

The cone we most frequently encounter is R ” , the nonnegative orthant, where 
R ” = K(I) with I being the identity matrix in R n * n . 

The cone K induces an order relation on R n by x^yoy-xeK, and 

x<yc-y-xeK°. The dual, or adjoint, cone is K* = (<f> e R”K<t>, x) > 0 V x € K}, 

which clearly satisfies the properties of a cone, and for a proper, simplicial cone, 

(K*)* = K If we define Ko as K\{0}, then x e K° (<j>, x) > 0 V <j> e K’ , and 

x e <3K <=> (<j), x) = 0 for some <j> e K * , where <9K is the boundary of K. 
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A continuous function f: D c: R n -» R n is quasimonotone nondecreasing in x 



relative to the cone K if for x,yeD,y-xe6K implies there exists a 4> e K * such 
that (<(>, y-x} = 0 and (<|), f(y) - fix)) > 0. This definition, from Eisner (1974), is fairly 
standard; however, for nonlinear functions some authors require x, y e K (see, for 
example, Heikkila, 1983). For a linear function fix) = Ax, A e R nxn , the quasimonotone 
nondecreasing property reduces to: x e 5K implies there exists a <j> e K* such that 
(<(), x) = 0 and (<(), Ax) > 0. 

An equivalent definition is f: R n — » R” is quasimonotone nondecreasing if and only 
if fj(x) is nondecreasing in xj for all i, j = 1 , . . ., n, i * j, so the following lemma 
characterizes quasimonotone nondecreasing linear operators. This result appears as an 
example in Lakshmikantham and Leela (1977b). 

Lemma A.l. A linear operator A e R n * n is quasimonotone nondecreasing 
relative to the nonnegative orthant if ay > 0 for all i * j. 

Proof. We prove this for n = 2. A similar argument works for n > 2. Let 
K - R^. — {(xi, x 2 )|xi > 0, x 2 > 0}. Then K* = K, since (<)), x) > 0 implies 
<> 1 X 1 + <j> 2 x 2 >0 for all xi, x 2 > 0, which in turn implies <J>i, 4>2 ^ 0. If f(x) = Ax, where 



A = 



an a 12 



a 21 a 22 



is quasimonotone nondecreasing relative to K, then for Xi = 0 or 



x 2 = 0, there exists a <J> e K * such that (<|), x) = 0 and (<]), Ax) > 0. Now 

(<t>, x) = <(>iXi + (J) 2 x 2 = 0 and (<(), Ax) = <j)i(anXi + a n x 2 ) + <t> 2 (a 2 iXi + a 22 x 2 ) > 0. If xi = 0, 

x 2 * 0, then <j) 2 x 2 = 0 so <J> 2 = 0. This imphes, since <()i > 0, that <t>iai 2 x 2 > 0, so ai 2 > 0. 
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Similarly, a 2 i > 0, so the linear functions which are quasimonotone nondecreasing relative 
to K are precisely those matrices with nonnegative off-diagonal entries. □ 

We use the following well-known results. 

Lemma A.2. Given a nonsingular matrix BeR" and the cone K(B), then 

K((B"') t ) = K(B)*. 

n n 

Proof. Let x 6 K((B -1 ) T ) and y e K(B) so that x = ^ ctj<|>i and y = ^ Phi, 

i=l i=l 

with a i3 pi > 0. Then <x, y> = ^ 2 Syaipj = > 0, so K((B _1 ) T ) c K(B)*. 

i=l j=l i=l 

n 

Similarly, let x g K((B _1 ) T ). Since the <f>i are independent, x = ^ yj<J>„ but there 

i=l 

exists at least one y, < 0. For y = bj e K(B), it follows that (x, y) = yj < 0, so x g K(B)*. 
Hence, K(B)* c K((B -1 ) T ), proving the le mma . □ 

Lemma A.3. Given a nonsingular matrix BeR" and the cone K(B), then 
K(B) cR" if and only if K(B _1 ) dR”, 

Proof. Let K(B)cR“ and let a e R” be arbitrary. Then Ba = xeR“ since 
x is a positive linear combination of the generators of K(B). Because B is invertible, it 
follows that a = B _1 x. As a e R ” is arbitrary and expressible as a nonnegative linear 
combination of the columns of B” 1 , we conclude that K(B -1 ) dR”. 

Conversely, assume K(B _1 ) dR” and let a e R” . Then there exists an 
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x e R ° such that B~'x = a, so x = Ba. Then B is nonnegative, otherwise it has some 
column bj which is not nonnegative, and letting a = [0 X , . . ., Oj.i, lj, Oj+i, . . ., 0 n ] T shows 
xgR°, which is a contradiction. This proves the lemma. □ 

Since K(B)cR" implies B is a nonnegative matrix, and since B T is also 
nonnegative, then K((B~') T ) 3 R ° as well. This gives a nice characterization of the 
inverses of nonnegative matrices. 

When we say that a matrix A e R nxn (or a linear operator A: R n — » R n ) is 
quasimonotone nondecreasing with respect to a cone K(B), we mean that C = B _1 AB is 
quasimonotone nondecreasing (with respect to R") even though A may not be. In 
Chapter HI we show this result as a theorem of Heikkila. 

Since quasimonotonicity is related to nonnegative matrices, we next present some 
background on nonnegative matrices. 

B. NONNEGATIVE MATRICES 

The study of nonnegative matrices began with Perron’s (1907) presentation of a 
theory of positive matrices. Frobenius (1908, 1909, 1912) immediately extended this 
theory to nonnegative matrices, and the general theory has carried their names ever since. 
Here we present only their basic results required for our application. 

A positive matrix A e R n>n has a^ > 0 for i, j = 1, . . ., n. A nonnegative matrix 
has a;j > 0. A matrix A is essentially nonnegative ( positive ) if and only if ajj > 0 (> 0) 
for all i^j. The following lemma justifies the term “essentially” nonnegative. 
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Lemma B.l. If a matrix A e R nxn is essentially nonnegative, then there exists a 
number r > 0 such that A = A + rl is nonnegative. Furthermore, the spectrum 
a( A ) = a(A) + r, and the eigenvectors of A and A are equal. 

Proof. The first statement is obvious. Let Ax = Xx. Then A x = (A + rl)x = 

Ax + rlx = Xx + rx = (X + r)x, completing the proof. □ 

Hence, we can shift the diagonal (and the spectrum) of an essentially nonnegative 
(positive) matrix to make it nonnegative (positive). 

A matrix A is reducible if it is permutation similar to a matrix 



Ap 



An A, 2 



0 A 



22 _ 



where An and A 2 2 are square (and of positive dimension). 



Otherwise, A is irreducible. A positive matrix is a special case of a nonnegative, 
irreducible matrix. The following theorem summarizes the basic results we need from the 
theory of nonnegative matrices. 

Theorem B.2. (Perron, 1907; Frobenius, 1912). A positive matrix A>0 has an 
eigenvalue equal to the spectral radius of A. Associated with this eigenvalue is a positive 
eigenvector. Furthermore, A has no other nonnegative eigenvectors. The same is true 
for a nonnegative, irreducible matrix A > 0. 

A nonnegative matrix A > 0 has a (real) eigenvalue equal to the spectral radius of 
A. Associated with this eigenvalue is a nonnegative eigenvector. □ 

For a proof of the above results see, for example, Gantmacher (1959), Varga 
(1962), or Horn and Johnson (1991). Since the transpose of a matrix with any of the 
above properties retains that particular property, the same conclusion can be made about 
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the left eigenvectors of A. Note that when comparing two matrices as in A > B, or two 
vectors as in x > y, we refer to component- wise majorization. 

From Lemma B. 1 and Theorem B.2 it follows that an essentially positive matrix 
has a real eigenvalue with greatest real part, with positive left and right eigenvectors 
associated with this eigenvalue, and no other nonnegative eigenvectors. Similar 
conclusions follow for essentially nonnegative and irreducible essentially nonnegative 
matrices. 

If a linear operator A: R n — > R" is quasimonotone nondecreasing, then the matrix 
A e R nxn is essentially nonnegative. Therefore, if we ask when a square, linear operator 
A on R" is quasimonotone nondecreasing with respect to a cone in R “ , we are asking 
the equivalent question of when the matrix A e R nxn is essentially nonnegative under a 
nonnegative change of basis. The next chapter introduces the motivating application, the 
theory of cone-valued vector Lyapunov functions. 
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III. CONE- VALUED VECTOR LYAPUNOV FUNCTIONS 



The motivating problem for determining the quasimonotonicity of a square matrix 
with respect to a nonnegative cone in R n comes from stability theory in differential 
equations. Specifically, we seek such a cone when using vector Lyapunov functions to 
determine the stability of an equilibrium in a dynamical system. This chapter presents this 
problem for linear comparison systems. 

We begin with a brief introduction to the theory of differential equations and 
stability. In order to develop the technique of vector Lyapunov functions, we also require 
results from the theory of differential inequalities. We present only the results from these 
fields needed to prove the major theorems of vector Lyapunov functions. In the case of 
well-known results, or those which will not be used directly, we present the results 
without proof 

A. DIFFERENTIAL EQUATIONS , INEQUALITIES, AND STABILITY 

We consider a system of first-order differential equations 

x' = fi;t, x), x(t 0 ) = xo, (A. 1) 

where x' = dx/dt, x e R m , f: I x D — > R m where I cz R and D cz R m are open sets 
containing t 0 and xo respectively, and f has certain continuity requirements. A 
solution x(t) to equation A. 1 is a differentiable function x: J cz R — » R m satisfying the 
differential equation dx/dt = f(t, x) and the initial condition x(t 0 ) = Xo. We also consider 
the autonomous system 

x' = fix), x(t 0 ) = xo, (A. 2) 
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where the vector field is independent of time. 

Picard (1890) and Lindelof (1894) showed that if f(t, x) is continuous on a closed 
set containing (t 0 , xo), and is uniformly Lipschitz continuous with respect to x on this 
set, then equation A. 1 has a unique solution x(t) on some set [t 0 , t 0 + a] for a > 0. 
Peano (1890) showed that without the assumption of Lipschitz continuity, then equation 
A.1 has at least one solution x(t) on t e [to, to + a]. Furthermore, a solution on an open 
subset of R m+1 containing (t 0 , xo) can be extended to the boundary of that set. Proofs of 
these results can be found, for example, in Hartman (1982), and these existence and 
uniqueness results are the basis for the theory of differential equations. 

To prove stability results about equilibria in equation A. 1, we also need results 
from the theory of differential inequalities. We extend the concept of the standard 
derivative used in equation A. 1 to the Dini derivative, defined as 

D_x(t) = liminf (l/h)[x(t + h) - x(t)]. This is one of four Dini derivatives, and we define 

h — >0 

the others as needed. We begin with scalar inequalities (m = 1 in equation A. 1). The 
following results are found, for example, in Lakshmikantham and Leela (1969a) or Walter 
(1970). 

Theorem A.l. Let E be an open (t, x)-set in R 2 and f e C[E, R]. Assume 

that v, w e C[[t 0 , t 0 + a), R] for some a > 0, and (t, v(t)), (t, w(t)) e E for 

t e [t 0 , t 0 + a). If v(t 0 ) < w (t 0 ) and if D_v(t) < f(t, v(t)) and D_w(t) > f(t, w(t)) for 

t e (t 0 , t 0 + a), then for t e [to, t 0 + a), v(t) < w(t). □ 

For the scalar differential equation A.1 (m = 1), we define the maximal solution on 
[t 0 , to + a) as the unique r(t) such that for every solution x(t) defined on [t 0 , t 0 + a), 
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x(t) < r(t) for t e [t 0 , t 0 + a). The minimal solution is defined analogously, and the 
following results presented for maximal solutions are also valid for minimal solutions, with 
appropriate changes. Under the hypothesis of Peano’s theorem, a maximal solution exists 
on the interval [to, t 0 + a], and as before, if one exists on an open set, it can also be 
extended to the boundary of that set. 

It is useful when employing vector Lyapunov functions to compare a solution to a 
maximal solution over an interval. The following comparison theorem is basic to this idea. 

Theorem A.2. Let f e C[E, R], where E is an open (t, x)-set in R 2 , and let 
[t 0 , t 0 + a) be the largest interval in which the maximal solution r(t) of equation A. 1 
exists. Let p e C[(to, to+ ct), R], (t, p(t)) e E for t 6 [t 0 , to + a), p(t 0 ) < xo, and for a 
fixed Dini derivative, Dp(t) < f(t, p(t)) on t e [t 0 , t 0 + a) except on possibly a countable 
subset of this interval. Then for t e [t 0 , t 0 + a), p(t) < r(t). □ 

The question of the existence of a maximal solution for a system of differential 
inequalities was solved by Wazewski (1950). He extended the above result to the case 
m > 1 in equation A. 1, and showed that a sufficient condition for a maximal solution to 
exist to equation A. 1 is that f(x) be non decreasing in xj for each t e [t 0 , to + a], a 
property which we defined in Chapter 13 as quasimonotone nondecreasing. This property 
was first recognized as being important to differential inequalities by Muller (1926) (see 
Walter, 1970). There are analogous definitions for increasing, decreasing, etc., as well as 
for mixed quasimonotone properties. 

Burton and Whybum (1952) used the mixed quasimonotone property to prove the 
existence of what are known as minimax solutions to differential equations. These are 
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solutions which, for example, are minimal in the first k components and maximal in the 
last m - k components. The existence of minimax solutions leads to a family of 
comparison theorems for differential inequalities which can be used to prove results about 
stability of equilibria using vector Lyapunov functions. These include results about 
stability, instability, and conditional stability. We present only the comparison theorem for 
stability, since we use stability for the motivating example of cone- valued vector 
Lyapunov functions. 

Theorem A.3. Let E be an open (t, x)-set in R m+1 , with f e C[E, R m ], 

Suppose that f is quasimonotone nondecreasing in x, and let [t 0 , to + a) be the largest 
interval of existence of the maximal solution r(t) of equation A. 1 . Let 
p e C[[t 0 , t 0 + a), R m ], (t, p(t)) e E for t e [t 0 , to + a), and on this interval, for the Dini 
derivative, D_p(t) < f(t, p(t)). Then p(to) < Xo implies p(t) < r(t) on this interval. □ 

This theorem applies for any Dini Derivative. The proof of this result, along with 
all the above results on differential inequalities, are in Lakshmikantham and Leela (1969a). 

We now present some basic definitions and concepts from stability theory which 
allow us to develop the technique of vector Lyapunov functions in its most simple setting. 

An equilibrium solution to equation A. 1 is a point E, such that f(t, £) = 0 for all 
t. Equation A. 1 admits the trivial solution E, = 0 if f(t, 0) = 0 for all t, and we frequently 
use this as our equilibrium as we can generally send an equilibrium point to the trivial 
solution via a change of coordinates. 

We wish to determine the stability of the equilibria of equation A. 1, and we begin 
with some basic definitions of stability. The trivial solution of equation A. 1 is stable if for 
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every neighborhood Ui of the origin and every to > 0, there is a neighborhood U 2 of the 
origin such that xo e U 2 implies x(t) e Ui for all t > to, where x(t) is a solution 
satisfying x(to) = xo. The trivial solution is asymptotically stable if it is stable and there 
exists a neighborhood U3 of the origin such that Xo e U3 implies that 
x(t) — > 0 as t — > 00 . 

We may express these definitions in terms of norms as follows. The trivial solution 
is stable if the solution x(t) with x(to) = xo is such that for every s > 0, to > 0, there 
exists a 5 > 0 such that ||xo|| < 5 implies ||x(t)|| < s for all t > to. Similarly, the trivial 
solution is asymptotically stable if there exists a 5 > 0 such that ||xo|| <5 implies that 
x(t) — » 0 as t —> 00 . 

The trivial solution in a linear system (when f(t, x) in equation A. 1 is a linear 
function of x, i.e., f(t, x) = A,(x) where A e R m<m ) is asymptotically stable when the 
eigenvalues of the matrix A all have negative real part for t > to- In this case we refer to 
the origin as a sink. In a nonlinear system, if the eigenvalues of the linearized system in a 
neighborhood of the equilibrium all have negative real part, then that equilibrium is a sink. 

If the origin is a sink in a linear system, then for all solutions x(t), ||x(t)|| is a 
strictly decreasing function of t, where ||*|| is the Euclidean norm. Similarly, it is a well- 
known result (see, e.g., Hirsch and Smale, 1974; or Arnold, 1989) that when the origin is 
asymptotically stable in a nonlinear system, there exists some norm such that ||x(t)|| is a 
decreasing function of t for all solutions x(t) starting sufficiently near the origin. 

These basic definitions of stability have numerous refinements. The two major 
classes of stability definitions are Lyapunov and Poisson, and among the logical statements 
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allowed in the definitions are 17,017,969 possible types of stability (see Bushaw, 1969). 
Many of these conditions are meaningless, and Massera (1949, 1956) discusses those 
which are most often studied in terms of Lyapunov functions. Habets and Peiffer (1973) 
extend these to a classification which encompasses all possible types of stability. 

The above definitions of stability can be generalized for nonautonomous systems 
(for example, when 5 depends on to) and they can be extended to conditional stability 
(instability) definitions for stable (unstable) manifolds through the point (to, 0). There is a 
well- developed theory of vector Lyapunov functions for conditional stability which uses 
the minima x solutions of Burton and Whybum, but we discuss only the most basic case. 
Furthermore, these stability results are local, and we do not present the conditions under 
which they are global. 

The definitions of stability and asymptotic stability used above rely on the norm of 
a solution being either bounded in any neighborhood of the trivial solution, or decreasing 
uniformly as t increases. Particularly for nonlinear systems, these conditions may be 
difficult to verify. Lyapunov recognized this, and in his 1892 dissertation he suggested 
that a function other than a norm could be used to determine stability (Lyapunov, 1907). 
Of course we now refer to such functions as Lyapunov functions, and while they are 
frequently difficult to find, they offer another powerful tool for classifying equilibria. 

Before introducing Lyapunov functions, we define the derivative of a function 
along a trajectory, in this case where the trajectory is a solution to the initial value 
problem For a solution x(t) to equation A. 1 and a real-valued function V(t, x), we 
define (using the terminology of Wiggins, 1996) the orbital derivative of V along x as 
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V (t, x) = limsup (l/h)[V(t + h, x + hf(t, x)) - V(t, x)]. We now present a version of 

h— >0* 

Lyapunov’s theorem from Hirsch and Smale (1974) for autonomous systems. A function 

V satisfying the hypothesis of the theorem is called a Lyapunov function. 

Theorem A.4. (Lyapunov, 1907). Let the autonomous system (equation A. 2) 
x'(t) = f(x) admit the trivial solution f(0) = 0. Let V(x) be a continuous function defined 
on a neighborhood U of the origin such that 0 e UcD. Let V be differentiable on 
U \{0}, with 

V(0) = 0, and V(x) >0 if x * 0. (A.3) 

Then 

V <0 on U\{0} (A. 4) 

implies the origin is stable. Furthermore, 

V < 0 on U\{0} (A. 5) 

implies the origin is asymptotically stable. 

Proof. Let U be any neighborhood of the origin, and choose 5 such that 
Bs(0) c U, where B s (0) is a 8-ball about the origin. Let a be the minimum value of V 
on the boundary of Bg(0), with a > 0 by equation A.3. Let Ui = {x e Bg(0)|V(x) < a}. 
Then no solution starting in Ui can meet the boundary of B s (0) since V is 
nonincreasing by equation A. 4. Therefore, the origin is stable. 

If equation A.5 holds as well then V is strictly decreasing on orbits in U\{0}. 

Let x(t) start in U]\{0} and suppose x(t n ) — > zo e B s (0) for t n — > oo. (By the 
compactness of Bg(0) such a sequence exists.) We now show that Zo = 0. By the 



19 



continuity of V and equation A. 5, V(x(t n )) -» V(zo) and V(x(t)) > V(zo) for all t > 0. If 
zo * 0, let z(t) be a solution starting at Zo. For any s > 0, V(z(s)) < V(zo), so for y(s) 
starting sufficiently near zo, V(y(s)) < V(zo). If we let y(0) = x(t n ) for sufficiently large 
n, then V(x(t„ + s)) < V(zo), which is a contradiction. Therefore, zo = 0 is the only 
possible limit point of {x(t)|t > 0}, and such a limit exists by the compactness of B§(0). 
This proves the theorem. □ 

While Lyapunov functions provide a useful technique for classifying equilibria, 
there is no known algorithm for finding them. In mechanical systems or electrical circuits, 
potential or total energy is frequently a candidate, but in general no technique works 
uniformly. Even after 100 years, finding a Lyapunov function is still more of an art than a 
science. The following example is from Hirsch and Smale (1974). 

Example A.5. Consider the system 
x'(t) = 2yz - 2y 



y'(t) = -xz + x 

z'(t) = -z 3 . 

The origin is clearly an equilibrium, but if we linearize the system near the origin via the 



Jacobian matrix then Jj 0 = 



0-2 0 
1 0 0 

0 0 0 



which has eigenvalues 0 and ±J2 i. We are 



unable to classify the origin via linearization, so we seek a Lyapunov function of the form 
V(x, y, z) = ax 2 + by 2 + cz 2 . Then 
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• • 

V = 2(axx' + byy' + czz') = 2[2axy(z - 1) - bxy(z -1) - cz 4 ]. Since we want V < 0, 

letting c = 1 and 2a = b gives V = -z 4 < 0, so V = x 2 + ly 2 + z~ is a Lyapunov 
function, and the origin is stable. However, the Lyapunov function is not strict since 

V = 0 on the xy-plane, and we cannot determine (with this choice of V) if the origin is 
asymptotically stable. □ 

While no known algorithm exists to find Lyapunov functions, they have proven 
very useful over the past one hundred years in determining the stability of equilibria. See, 
for example, Arnold (1989), Bailey (1966), Bhatia and Szego (1970), La Salle and 
Lefschetz (1961), Lehnigk (1966), and Yoshizawa (1966). 

In 1962, Bellman recognized that the requirement that a Lyapunov function be 
real- valued was too restrictive, and he proposed that a vector-valued function might 
provide more flexibility. He presented his results for square, linear systems and special 
cases of nonlinear systems. The following lemma is a special case of Theorem A. 3 for 
linear systems, but we present it in its entirety as it provides necessary and sufficient 
conditions for the solution of a differential equation to be majorized by a solution of a 
differential inequality. In this le mma . A = [aij] e R mxm . The le mma is from Beckenbach 
and Be llman (1965). 

Lemma A.6. Let the system of differential equations 

dx/dt = Ax, x(0) = xo 

have a solution x(t). Let y(t) satisfy the differential inequality dy/dt > Ay, y(0) = Xo. 
Then y(t)>x(t) for t>0 if and only if a;j > 0 for i^j. 
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Proof. The system dy/dt = Ax + g(t), y(0) = xo has a solution of the form 



y(t) = xoe At + J" e A(l ' S) g(s) ds, so we need to know when the elements of e At are 

0 

nonnegative for t > 0. Since e At = I + At + A 2 t 2 /2! + A 3 t 3 /3! + . . ., for small positive t 
e At « I + At so that y(t) > x(t) imphes a;j > 0. 

Conversely, assume ay > 0 (the case ay > 0 follows via a limiting procedure) and 
let e At = ( e At/N ) N . For a fixed t, e At/N = I + A(t/N) + A 2 (t/N) 2 /2! + . . . is positive (in the 
sense that all elements are positive for N sufficiently large). Since the product of positive 
matrices is positive, e At is positive if ay > 0 for i^j. This proves the lemma. □ 

The condition that ay > 0 for i^j is precisely the condition of f(x) = Ax being 
quasimonotone nondecreasing in Theorem A. 3 for a linear function. Also, since a linear 
function is uniformly Lipschitz continuous, Picard’s theorem applies and the maximal and 
minimal solution are the unique solution x(t) = xoe At . 

We now state a stability theorem for vector Lyapunov functions, where we 
dete rmin e the stability of the trivial equilibrium solution to equation A. 1. 

Theorem A.7. (Lakshmikantham, 1965). Let g e C[R+ x R” , R n ], g(t, 0) = 0 
for all t, let g(t, u) be quasimonotone nondecreasing in u for each t e R+, and let 

u ' (t) = g(t, u), u(t 0 ) = Uo > 0. (A.6) 

n 

Let V e C[R+x D, R" ] and V(t, x) be locally Lipschitz in x, with T] V,(t, x) — » 0 as 

i=l 

||x|| — » 0 for each t; and let f e C[R+ x D, R m ], f(t, 0) = 0, and for (t, x) € R+ x D, 

V (t, x) < g(t, V(t, x)). 
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Then the stability of the trivial solution u = 0 of equation A.6 implies the stability of the 
trivial solution x = 0 of equation A. 1 . □ 

In the more general theorem of Lakshmikantham (1965), various types of stability 
and conditional stability are proven as well. Instead of presenting the proof here, we 
present instead an example from Lakshmikantham (1974) that demonstrates the usefulness 
of a vector Lyapunov function over that of a scalar Lyapunov function. 

Example A.8. Consider the system 

dx/dt = e”‘x + ysint - (x 3 + xy 2 )sin 2 t, (A. 7a) 

dy/dt = xsint + e~‘y - (x 2 y + y 3 )sin 2 t. (A. 7b) 

We would like to determine the stability of the trivial solution to this system, so we 
attempt to find a scalar Lyapunov function, trying V(x, y) = x 2 + y 2 . The best bound we 



can achieve with this choice of V is V < 2(e -t + |sint|)V. Clearly the trivial solution of 
du/dt = 2(e _t + |sint|)u is not stable, so we cannot conclude anything about the stability of 



the trivial solution of equation A. 7 with this particular choice of scalar Lyapunov function. 



However, if we choose V(x, y) = 



V](x, y) 
v 2 (x,y). 



= (1/2) 



(x + y) 2 

(x-y ) 2 



as a vector 



Lyapunov function, then V i < 2(e _t + sint)Vi, and V 2 < 2(e _t - sint)V 2 . Since the trivial 
solution u = 0, w = 0 to the system 

du/dt = 2(e _t + sint)u 
dw/dt = 2(e _t - sint)w 

is clearly stable, the trivial solution x = 0, y = 0 to the system A. 7 is also stable. □ 
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While we used the quantity ^ V;(t, x) as a measure, we could have used some 

i=l 

other measure such as max Vj(t, x) or Q(Vi(t, x), . . V n (t, x)) where Q: R R+ 

i 

and Q(u) is monotone nondecreasing in u. 

Vector Lyapunov functions offer a lot of flexibility when they can be found, but 
their true versatility is apparent in the case of conditional stability; i.e., when the stable 
manifold of the equilibrium has dimension less than m. Here we present the most basic 
results of this theory in terms of conditional stability, which are from Lakshmikantham 
(1965). For completeness we review some of the terminology. 

We wish to determine the stability of the (generally trivial) equilibrium solution to 
equation A. 1, 

x'(t) = f(t, x), x(t 0 ) = xo, t 0 > 0, (A. 1 ) 

where x 6 R m and f(t, x) is defined and continuous in R+ x R m . The trivial solution 
x = 0, where f(t, 0)s0 for all t > 0, is conditionally equistable (where 8 in our 
definition of stability depends on s and t 0 ) if and only if there exists a manifold M (m _k) 
through the origin of dimension (m - k) such that for each s > 0 and each t 0 > 0 there 
exists a positive function 8 = 8(to, e) that is continuous in t 0 for each s, where ||x(t 0 )[| ^ 
8 and x(to) e M( m _k) imply ||x(t)|| < s for all t > to. 

We consider the comparison system 

r'(t) = g(t, r), r(t 0 ) = r 0 > 0 (A. 8) 
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where g(t, r) is quasimonotone nondecreasing in r e R n . When g has this property, the 
solution r(t; to, r 0 ) is maximal in the sense of component- wise majorization. We consider 
solutions of the form 

r i(t; to, ro) = 0 (i = 1, 2, . . ., k) (A.9a) 

ri(t; t 0 , r 0 ) >0 (i = k+1, . . n) (A. 9b) 

where m - k corresponds to the dimension of the manifold M (m _ky 

The comparison system A. 8 is equistable (where 5 in our definition of stability 
depends on s and to) if and only if for all z > 0, t 0 > 0, there exists a positive function 

n 

5 = 8(to, s), continuous in to for each s, such that ^ r 0 <5 implies 

i=k+l 

n 

^ r;(t; t 0 , r 0 ) < s for t > to. Other stabihty and conditional stability properties are 

i=i 

formulated similarly. 

The vector Lyapunov function V(t, x) e C[R+ x R m , R " ] is locally Lipschitz 
continuous in x, and its orbital derivative along the trajectories x(t) of equation A. 1 
satisfies 

V(t,x)<g(t,V(t,x)) (A. 10) 

where x(t) is any solution such that V(to, xo) < ro. By Theorem A. 3 this ensures that 
V(t, x(t)) < r(t; t 0 , r 0 ) for t > t 0 . 

The stable manifold of the origin, M^-k), is defined by the set of points for which 

Vi(t, x) s 0 for i = 1, 2, . . ., k < m (A. 1 1) 
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The bound for the measure, b(r) is a continuous function, nondecreasing in r, 
such that 

b(r) > 0 for r > 0, and b(||x||) < Vi(t, x). (A. 12) 

i=l 

We might also have the bound satisfy b(r) -» oo as r — » oo, which is required for certain 
types of stability. 

While we have used the component- wise sum of the Lyapunov function as the 
measure, there are several equivalent choices we could have selected. In terms of this 
measure, we may require, for example, 

n 

^ V;(t, x) — > 0 as ||x|| — > 0 for each t > 0. (A. 13) 

i=l 

We now state the theorem on vector Lyapunov functions and conditional stability. 
Theorem A.9. (Lakshmikantham, 1965). Let assumptions A. 10 through A. 13 
hold. If the solution A. 9 to equation A. 8 is conditionally equistable, then so is the 
equilibrium solution of equation A. 1 

Proof. Let s > 0. If ||x|| = e, then from assumption A. 12, we have 

b<s)S £ Vj(t,x). (A. 14) 

i=l 

If the stability property holds, given b(s) > 0 and to > 0, there exists a positive function 
8 = 8(to, s), continuous in to for each s, such that 

D 

X r ; (t; t 0 , r 0 ) < b(s) for t > t 0 , (A. 15) 

i=l 
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n 



provided ^ r 0 <8. Let x(t) be any solution of equation A. 1 . 

i=k+l 


Then it follows from 


assumption A. 10 that 




n n 

2 Vj(t, x(t)) < 2 rj(t; t 0 , r 0 ) for t > t 0 

1=1 i=l 


(A. 16) 


whenever 




S V<to, x(to)) < 2 V 

i=l i=l 


(A. 17) 


Now choose r 0 (i = 1, 2, . . n) to satisfy 




r 0 =0 for i = 1, 2, . . k and 


(A 18) 


£ r 0j -8. 


(A 19) 



i=k+l 



By equation A. 17 and since V;(t, x) > 0 for i = 1, 2, . . n, equation A. 1 8 implies that 



x(t 0 ) e M( m _ k) because of assumption A. 1 1. From the monotonic property of b(r), 
assumption A. 12, and equations A. 17 and A. 19, we deduce that ||x(to)|| < b'^S) = Si. 



By assumption A. 13, there exists a 82 = 8 2 (t 0 , 8) such that sup ^ V,(t 0 , x(to)) < 8. 

Bx«o)l|sS ; i=i 



Let 83 = min{8i, 8 2 }. It then follows from the choices of r 0 and S3 that x(t 0 ) e M< n . k ) 



and ||x(t 0 )|| < 83 implies every solution x(t) satisfies equation A. 16. Suppose, if possible, 
that a solution x(t) of equation A. 1 satisfying that x(t 0 ) e M( n .k) and ]]x(to)]| ^ 83 is such 
that ||x(x)|| = s for some x > t 0 . Then from assumption A. 12 and equations A. 14, A. 15, 



and A. 16 follows the contradiction b(s) < 2 Vi(x, x(x))<^ r;(x; t 0 , r 0 ) < b(s), 

i=l i=l 
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which proves the theorem. □ 

Vector Lyapunov functions can also be used to determine certain boundedness 
properties of solutions of equation A. 1, and boundedness theorems analogous to the 
above theorems for stability are given in Lakshmikantham (1965). We now give an 
example from Lakshmikantham (1965) of the usefulness of a vector Lyapunov function in 
determining the conditional stability of an equihbrim 

Example A.10. Consider the system of differential equations 

dx/dt = (1 + cos(t))x + (1 - cos(t))y + (cos(t) - l)z 

dy/dt = (1 - e“')x + ( 1 + e~‘)y + (e"‘ - l)z 

dz/dt = (cos(t) - e‘‘)x + (e -t - cos(t))y + (e"‘ + cos(t))z. 

We select as the Lyapunov function 





Vj(t,x,y,z) 




(x + y-z) 2 


V(t, X, y, z) = 


V 2 (t,x,y,z) 


= 


(x-y + z) 2 




_Vj(t,x,y,z)_ 




(-x + y + z) 2 



3 

Since V; = x 2 + y 2 + z 2 + (x - y) 2 + (y - z) 2 + (z - x) 2 , we may use as the positive, 

i=l 

monotone function b(x, y, z) the square of the Euclidean norm, since 

3 

b(x, y, z) = x 2 + y 2 + z 2 < ^ Vj(t, x, y, z), so condition A. 12 is satisfied. Since 

i=l 

V i(t, x, y, z) = 2(x + y - z)( dx/dt + dy/dt - dz/dt) = 4(x + y - z) 2 

V 2 (t, x, y, z) = 2(x - y + zX dx/dt - dy/dt + dz/dt) = 4cos(t)(x - y + z) 2 

V 3 (t, x, y, z) = 2(-x + y + zX -dx/dt + dy/dt + dz/dt) = 4e~*(-x + y + z) 2 . 
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the comparison system is defined by 



gi(t, V,,V 2 , V 3 ) = 4Vi(t, x, y, z) 
g 2 (t, Vi, V 2 , V 3 ) = 4cos(t)V 2 (t, X, y, z) 
g 3 (t,V 1 ,V 2 ,V 3 ) = 4e- t V 3 (t > x,y,z). 
The comparison system is therefore 



r/(t) 

t 




Sl(t> r i ,h > r 3 ) 




4r, 


r 2> (0 


= r '(t) = g(t, r) = 


g 2 (t, r i , r 2 , r 3 ) 


= 


4 cos(t)r 2 


r 3 '(t) 




_g3(t,ri,r 2 ,r 3 )_ 




4e_ ‘ r 3 



Since each gj is only a fimction of r;, the fimction g is trivially quasimonotone 

nondecreasing, and satisfies v <g(t,V). 

Since the trivial solution to rf(t) = gl (t, r) is clearly unstable, we must choose 
k = 1, so the initial condition is r 0i = 0, and ri(t; to, 0) = 0. This forces the condition 

Vi(t, x, y, z) = 0, which defines the two-dimensional manif old M( 3 _i) by the condition 
(x + y - z) 2 = 0 to be the plane x + y = z. The solution to the comparison system is 
r i(t) = 0, 

r 2 (t) = r 0; exp[4(sin(t) - sin(t 0 ))]. 



r 3 (t) = r 0j expH^-e'* 0 )]. 

3 

If s > 0, then for 8 = min{s/(2e 8 ), s/(2exp(4e ~‘° )}, r 0j + r 0j <5 imphes ^ r;(t) < e 

i— 1 

for t > t 0 . Since the conditional equistability condition is satisfied, the trivial solution to 
the original system is equistable if the initial point is on the manifold x + y = z. □ 
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The technique of vector Lyapunov functions is a proven technique in stability 
analysis in many fields, particularly for large scale dynamical systems and interconnected 
systems. See, for example, Michel and Miller (1977), Siljak (1978), Lakshmikantham, 
Matrosov, and Sivasundaram (1991), and Vidyasagar (1993). 

B. CONE-VALUED VECTOR LYAPUNOV FUNCTIONS 

While vector Lyapunov functions provide a great amount of flexibility over scalar 
Lyapunov functions, Lakshmikantham (1974) noted that the quasimonotonicity of the 
comparison system is not a necessary condition for the system to be stable. In particular, 
for a linear system, a matrix can still be a stability matrix (have all eigenvalues with real 
part less than zero) without having all the off-diagonal elements nonnegative. 

By requiring the vector Lyapunov function to have each V; > 0, we have restricted 
it to the cone R “ . Lakshmikantham and Leela (1977a) investigated the possibility of 
selecting a cone other than R “ to overcome this limitation. The following theorems 
from Lakshmikantham and Leela (1977b) extend some of the previous results to cones 
other than R ” . The first is the comparison principle corresponding to Theorem A.3 
through the cone K. 

Theorem B.l. (Lakshmikantham and Leela, 1977b). Let f e C[R+ x R n , R n ] be 
quasimonotone nondecreasing in x relative to K for each t e R+, and let [t 0 , go), to > 0, 
be the largest interval of existence for the maximal solution r(t; to, Xo) of equation A. 1 
relative to K. Further, let p e C[R+, R n ] and D_p(t) < f(t, p(t)) for t > to. Then 
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p(t 0 ) | xo implies p(t) | r(t; t 0 , xo) for t > t 0 . □ 

The quasimonotonicity of f(t, x) in x relative to the cone P does not necessarily 
imply the quasimonotonicity of f(t, x) in x relative to the cone Q when 
P c Q. However, if P c Q, then the order relations relative to P do imply the same 
order relations relative to Q. We prove these observations in Section D, and from them 
comes the following corollary. 

Corollary B.2. Let P and Q be two cones in R n such that PcQ. Let the 
assumptions of Theorem B. 2 hold with K replaced by P. Then p(t 0 ) < Xo implies 

p(t) < r(t; t 0 , xo) for t > t 0 . □ 

Q 

We now state the comparison theorems and stability results for cone-valued vector 
Lyapunov functions from Lakshmikantham and Leela (1977b). We begin with the system 
of differential equations 

x'(t) = fi;t, x), x(t 0 ) = xo, (A. 1) 

where f e C[R+ x D, R m ]. If K is a cone in R", n < m, and the cone- valued vector 
Lyapunov function V € C[R+ x D, K], we define for (t, x) € R+ x D, the orbital Dini 
derivative as D + V(t, x) = limsup (l/h)[V(t + h, x + hf(t, x)) - V(t, x)]. The first result 

h— 

follows. 

Theorem B.3. (Lakshmikantham and Leela, 1977b). Assume that V(t, x) 
satisfies a local Lipschitz condition in x relative to K and for (t, x) e R+ x D, 

D*V(t, x) | g(t, V(t, x)), (B.l) 
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where g e C[R + x K, R n ], and g(t, u) is quasimonotone in u with respect to K for 
each t € R+. If r(t; t 0 , uo) is the maximal solution of u'(t) = g(t, u), u(0) = uo relative to 
K, and x(t; t 0 , Xo) is any solution of equation A. 1 such that V(t 0 , xo) < uo, then on the 

common interval of existence 

V(t, x(t; t 0 , xo)) | r(t; to, u 0 ). (B.2) 

Proof. Let x(t) = x(t; to, xo) be any solution as above. Set p(t) = V(t, x(t)). For 
small h > 0, since V(t, x) is locally Lipschitz in x with respect to K, then 
p(t + h) - p(t) | L||x(t + h) - x(t) - hf(t, x(t))|| + V(t + h, x(t) + hf(t, x(t))) - V(t, x(t)). 

From this and equation B. 1 follows the inequality D + p(t) ^ g(t, p(t)). Applying Theorem 

B. 1 gives the conclusion B.2. □ 

The following variant of this theorem offers more flexibility in applications, and its 
proof follows from Corollary B.2. 

Theorem B.4. (Lakshmikantham and Leela, 1977b). Let P and Q be two cones 
in R n such that P c Q. Suppose that V e C[R+ x D, Q], V(t, x) satisfies a local 
Lipschitz condition in x relative to P, and D + V(t, x) < g(t, V(t, x)) for (t, x) e R+ x D, 

P 

where g e C[R+ x D, R n ] and g(t, u) is quasimonotone nondecreasing in u relative to 
P for each t € R+. If r(t; t 0 , uo) is the maximal solution as in Theorem B.3 relative to P 
and x(t; t 0 , xo) is any solution of equation A. 1 such that V(t 0 , xo) < Uo, then 

P 

V(t, x(t; t 0 , xo)) < r(t; t 0 , Xo) (B.3) 

Q 
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on the common interval of existence of r(t; t 0 , uo) and x(t; to, uo). When Q = R " , 
relation B.3 implies the componentwise inequalities V(t, x(t; to, Xo)) < r(t; t 0 , uo)- □ 

We now state the first stability results for cone-valued vector Lyapunov functions. 
The first theorem extends the result of Theorem A. 7 for equistability. Its proof follows 
from that of Theorem A.9 and Theorem B.3. 

Theorem B.5. (Lakshmikantham and Leela, 1977b). Let the assumptions of 
Theorem B.3 hold. Let f(t, 0) = 0 and g(t, 0) = 0, and assume that for some 
<j>o e K; and (t, x) e R+ x D, 

b(||x||) < <<|>o, V(t, x)> < a(t, ||x||), (B.4) 

where a e C[R+ x [0, p), R+], b e C[[0, p), R+], a(t, 0) = 0, b(0) = 0, a(t, u) and b(u) are 
increasing in u, and p = sup{||x|||x e D}. Let the trivial solution u = 0 of u'(t) = g(t, u) 
be <j) 0 -equistable, that is, given s > 0, to > 0, there exists a 8 = 8(to, s) > 0 such that 
(<j)o, uo) < 8 imphes (<J) 0 , r(t; to, uo)) < e for t>t 0 . Then the trivial solution x = 0 of 
equation A. 1 is equistable. □ 

The following theorem again increases our flexibility when employing cone-valued 
vector Lyapunov functions. The proof follows from those of Theorems A.9 and B.4. 

Theorem B.6. (Lakshmikantham and Leela, 1977b). Let the assumptions of 
Theorem B.4 hold, and let f(t, 0) s 0 and g(t, 0) = 0. Assume equation B.4 is satisfied 
for some <j> 0 e Q ' 0 and that the trivial solution usOof u'(t) = g(t, u) is <j>o-equistable, 
with <J>o 6 Q j . Then the trivial solution x = 0 of equation A. 1 is equistable. □ 



If K = R“ , then <t> 0 = [1, 1, . . 1] T , and Theorem B.5 reduces to Theorem A.7. 



In this case, condition B.4 becomes b(||x||) < V,(t, x) < a(t, ||x||). We could again use 

i=l 

measures other than the component-wise sum of the vector Lyapunov function. In 
Theorem B.6, if P c Q = R ” , we can remove the requirement for quasimonotonicity if 

we select an appropriate cone P depending on the nature of g(t, u) as shown in the 
following example from Lakshmikantham and Leela (1977b). 

Example B.7. Consider the system 

ui'(t) = anUi + ai 2 u 2 = gi(t, ui, u 2 ), ui(t 0 ) = u 0 i, (B.5a) 

u 2 '(t) = a 2 iui + a 22 u 2 = g 2 (t, ui, u 2 ), u 2 (t 0 ) = uo 2 . (B. 5b) 

If Q = R 2 ,and a 2 i and ai 2 are not nonnegative, then the function g(t, u) is not 
quasimonotone nondecreasing in u = (ui, u 2 ) relative to Q. Hence, the differential 
inequalities 

D + V i(t, x) < gl (t, V,(t, x), V 2 (t, x)), (B.6a) 

D^t, x) < g 2 (t, V,(t, x), V 2 (t, x)) (B.6b) 

do not yield the componentwise estimates of V(t, x(t)) in terms of the solution of 
equation B.5. However, if there exist two numbers a and P such that 0 < P < a and 

a 2 a 2 i + aa 22 > aau + ai 2 , (B.7a) 

P 2 a 2 i + Pa 22 > Pan + a J2 , (B.7b) 

then these conditions hold with no restriction on the nonnegativity of a 2 i and ai 2 . Let 
the cone PcQ = R 2 + be defined by P = (u e R \ |Pu 2 < Ui < au 2 }. The boundaries for 
this cone are au 2 = ui and Pu 2 = ui. On the first boxmdary, let (J) = (-l/a, 1) so 
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<(-l/a, 1), (ui, ui/a)> = 0, and <(-l/a, 1), (a u ui + a^ui/a, a 2 iUi + a 2 2 U]/a)) > 0 for all 
u * 0. This reduces to condition B. 7a, and we similarly obtain condition B. 7b. Thus, if 
the inequalities B.6 are relative to P, we obtain component- wise estimates on V as 

V ; (t, x(t)) < r ; (t; t 0 , V(t 0 , xo)) (B. 8) 

by Theorem B.4. If a ]2 , a 2 i > 0, then equation B.8 is the one we would obtain through the 
standard method of vector Lyapunov functions. □ 

While it appears to be a nontrivial exercise to construct an appropriate cone for a 
system of differential equations, this method still carries much merit as it provides a further 
increase in our flexibility to determine the stability of equilibria. 

These are the most basic results from the theory of cone-valued vector Lyapunov 
functions we require to motivate the problem of finding a nonnegative cone with respect 
to which a linear differential operator on R n is quasimonotone nondecreasing. Further 
directions in the theoiy of vector Lyapunov functions have been explored recently. In 
particular, using higher derivatives of vector Lyapunov functions, Koksal and 
Lakshmikantham (1996) showed how to find a particular cone with respect to which a 
given comparison system (resulting from taking higher derivatives of a decrescent 
Lyapunov function) is quasimonotone. 

Furthermore, there has been much research on nonlinear comparison systems. See, 
for example, Hatvani (1984), Deimling and Lakshmikantham (1990), Lakshmikantham, 
Leela, and Ram Mohan Rao (1991), and Lakshmikantham and Papageorgiou (1994). 
However, we focus on finding cones for vector Lyapunov functions when the comparison 
system is linear. 
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C. THE QUASIMONOTONICITY OF LINEAR DIFFERENTIAL SYSTEMS 

The problem of finding a cone in R" with respect to which a given linear operator 
is quasimonotone nondecreasing (or essentially nonnegative) and many related problems 
have been addressed in various forms in the literature in recent years. For example. 

Van dergrafi (1968) following Birkoff (1967) gave sufficient conditions for such a cone to 
exist without further restrictions on the cone K. However, Heikkila first addressed this 
problem for the application of vector Lyapunov functions, with the requirement that the 
cone be proper, simplicial, and nonnegative. He stated this to be sufficient for the 
application of vector Lyapunov functions, and in the next section, we justify this claim. In 
this section, we summarize Heikkila’ s results. 

Consider a linear mapping A e R nxn and a cone K generated by n vectors b, in 
R n such that 

K = {wibi + . . . + w n b n |wj > 0}. (C. 1) 

For C. 1 to define a proper, simplicial cone in R", it is necessary and sufficient that the bi 
be linearly independent. The first result gives necessary and sufficient conditions for a 
linear operator to be quasimonotone with respect to a cone in R n . 

Theorem C.l. (Heikkila, 1980). A linear mapping A: R n — > R" is 
quasimonotone nondecreasing relative to the cone generated by a basis {bi, . . ., b n } of 
R n if and only if the matrix of A relative to this basis has all off-diagonal elements 
nonnegative. 

Proof. Let A = [a,j] and C = [cy] be the n x n matrices of A relative to the 
standard basis of R" and the basis (bi, . . ., b n } respectively. Then C = B -1 AB, where 
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B e R nxn has bj as its jth column. Let <]) J be the ith row of B \ so 

<4>i, bj) = Sij, (C.2) 

and 

<<j>i, Abj) = cij. (C.3) 

Equations C.2 and C.3 imply that (j)j e for each i = 1, . . ., n. 

n 

Assume cjj >0 for i * j, and let u e 3K be given. Then if u = ^ wjbj, at least 

j=i 

n 

one coefficient w, must be zero. Then <4>i, u> = ^ Wj(<j)j, bj) = 0 and 

j=i 

n n 

<<]);, Au)=^] Wj(<J>i, Abj) = ^ WjCjj>0 from equations C.2 and C.3. Thus, the 

j=i j=i 

condition Cjj > 0 for i ^ j impbes the quasimono tonicity of A relative to K. 

Conversely, let at least one off-diagonal element of C, say c,j, is negative. Then 
by equation C.3, ((j);, Abj) < 0. Since the mapping u — > (<J>;, Au) is continuous, there 
exists a 8 > 0 such that ||u - bjll < 5 impbes 

<<t>i, Au) < 0. (C.4) 

D 

In particular, equation C.4 holds for u e 5K given by u = ^ Wkbk, with w; = 0, 

k=l 

Wj = 1, and for k =£ i, j, Wk = 8/(n||bk||). Since all the coefficients w k , except w;, are 
positive, then for any b e K* for which (b, u> = 0, it follows that (b, bk> = 0 for aU 
k^i. Thus, b must be of the form b = a<j)i for a > 0, so by equation C.4, 



37 



(b, Au) = a{<J)j, Au) < 0. This shows A is not quasixnonotone nondecreasing relative to 
K unless Cjj >0 for i ^ j, completing the proof. □ 

Salzmann (1972) proved the same result for nonnegative matrices, and he 
presented it in terms of positive operators on simplicial cones. Furthermore, in his 
discussion of matrices with invariant cones, Vandergraft (1968) proved a s imilar result. An 
immediate consequence is the following. 

Corollary C.2. (Heikkila, 1980). If a linear mapping A: R n — » R n has n 
linearly independent eigenvectors bi, . . .b n , then A is quasimonotone nondecreasing 
relative to the cone generated by these eigenvectors. 

Proof. In this case, the matrix of A relative to {bi, . . ., b n } is diagonal. □ 

We can also express a necessary condition for monotonicity in terms of the 
eigenvectors of A. This result follows from the result from Perron ( 1907) and Frobenius 
(1912) on the theory of nonnegative matrices, which states that a nonnegative matrix has a 
nonnegative eigenvector corresponding to a real eigenvalue of greatest modulus (see 
Theorem n.B.2). 

Corollary C.3. (Heikkila, 1980). If a linear mapping A: R n — » R" is 
quasimonotone nondecreasing relative to a cone of type C. 1, then A has a non-zero 
eigenvector in this cone, and hence a real eigenvalue. 

Proof. If A<, = B -1 AB is quasimonotone nondecreasing where B = [bi, . . ., b n ], 
choose a positive number r so large that C = A<, + rl is nonnegative. Then C has a 
nonnegative eigenvector xi > 0, so Aq (and hence A) has an eigenvector 
yi = Bx! e K(B). □ 
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We now consider linear differential systems. By the quasimonotonicity of the 
linear system 

u'(t) = A(t)u, t>0, (C.5) 

relative to a cone K we mean that for each t > 0 the linear mapping A(t) is 
quasimonotone nondecreasing relative to K. The next result for cones contained in R “ 
follows from Corollaries C.2 and C.3. 

Corollary C.4. (Heikkila, 1980). For the quasimonotonicity of the system C.5 
relative to some cone in R " it is necessary that for each t > 0 at least one eigenvector 
of A(t) belongs to R“, and sufficient that A(t) has n linearly independent eigenvectors 
in R “ , which do not depend on t. □ 

Given a nonsingular matrix B e R nxn , the mapping u = Bv transforms equation 
C.5 into the equivalent form 

v'(t) = B _1 A(t)Bv, t>0. (C.6) 

But B _1 A(t)B is the matrix of the linear mapping A(t) in the basis of R n formed by the 
col umn vectors of B. Moreover, if B is nonnegative, these vectors belong to R“ . 

Hence, the next result follows from Theorem C. 1 . 

Theorem C.5. (Heikkila, 1980). There exists a cone K in R“ generated by n 
linearly independent vectors, relative to which the system C.5 is quasimonotone 
nondecreasing if and only if there exists a nonsingular matrix B such that the off-diagonal 
elements of the coefficient matrix of the system C.6 are nonnegative. Moreover, the 
col umn vectors of B generate K. □ 
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This special case of considering cones in R ” lets us apply Theorem B.6 when the 
vector Lyapunov function satisfies a differential equality. 

D. STABILITY THROUGH LINEAR COMPARISON SYSTEMS 

The results of Heikkila (1980) presented in Section C address the problem of 
finding a cone with respect to which a linear comparison system is quasimonotone, and 
provide a necessary condition and a sufficient, condition for, such a cone to exist. Koksal 
and Lakshmikantham (1996) constructed cones for comparison systems generated by 
taking derivatives of Lyapunov functions, and Koksal and Fausett (1995) extended these 
results using generalized eigenvectors. 

Heikkila (1980) states (without proof) that for stability results it is sufficient to 
consider nonnegative cones. This section gives our proof of this result, and shows it is 
sufficient to consider square, linear comparison systems. 

We wish to determine the stability properties of the origin as a solution to 

x'(t) = f(t, x) (A. 1) 

where f e C[R+ xR% R m ] and f(t, 0) = 0. To do this we seek a vector Lyapunov 
function V(t, x) € C[R+ x R m , K] which is locally Lipschitz in x with respect to a cone 
K c R n , and which satisfies the differential inequality 

V | g(t, V) (D.l) 

where V is the orbital derivative of V along trajectories x(t). If g(t, u) is 
quasimonotone nondecreasing with respect to K, the stability of the trivial solution of 
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U'(t) = g(t, U) 



(D.2) 



may be used to determine the stability of the origin x(t) = 0 in equation A. 1. We require 
certain conditions on the Lyapunov function V for this result to hold, for example, 

b(l|x||)<£ |VJ<a(||x||) (D.3) 

i-1 

where a and b are as in Theorem B. 5. 

We let the comparison system A.3 be a square, linear system, or g(t, u) = A(t)u 
where A: R+ x R n — » R”, or for all t, A e R n * n . Let V e R " so the quasimonotonicity 
of A(t) implies its off-diagonal elements are nonnegative, and the inequality D. 1 is 
component-wise majorization (otherwise we can change bases so V e R " ). 

If A(t) is a stability matrix, but does not possess the required quasimonotonicity 
properties, we seek a cone KcR" other than R " for which inequality D. 1 holds, 
relative to which A(t) is quasimonotone nondecreasing, and which contains V(t, x) 
while maintaining its locally Lipschitz property. The following lemma discusses the effect 
of a linear transformation on an inequality. 

Lemma D.l. If x < y, then for a nonsingular matrix B e R nxn , Bx < By, 

r; k(b> 

where K(B) <z R n is the cone generated by the columns of B. 

Proof. Since Bx < By e=> By - Bx e K(B) and (y - x) e R ", then from the 

K(B) 

□ 

linearity of B, B(y - x) e K(B) = (u = ^ a;bi|ai > 0, bj is a column of B}. Letting 

i=l 

oti = (y — x)j > 0 completes the proof. □ 



41 



We can transform a differential system v'(t) = Av, with a nonsingular matrix B, 
where y = Bv, so y'(t) = BAB”'y. Such a transformation is called a Lyapunov 
Transformation and the stability of the origin in one system implies the stability of the 
origin in the other. In this case, we say the systems are equivalent in the sense of 
Lyapunov (see Gantmacher, 1959). 

We now show that a simple coordinate transformation cannot be used to construct 
a cone which satisfies the hypothesis of Theorem B.6 when the original system does not 
possess the required quasimonotonicity property. 

Proposition D.2. Assume a vector Lyapunov function satisfies the inequality 

V < A(t)V, (D.4) 

r; 

but that A(t) is not quasimonotone nondecreasing with respect to R ” . Then there does 
not exist a Lyapunov transformation Y = B _1 V for which Y < C(t)Y and C(t) is 

K(B~' ) 

quasimonotone nondecreasing with respect to K(B _1 ). 

Proof. If such a transformation existed, then with C(t) = B -1 A(t)B, multiplying 
inequality D.4 by B _1 and applying Lemma D. 1 gives 

Y = B _1 V < B"'A(t)V = B -1 BC(t)B -1 V = C(t)Y. 

K(B-‘) 

However, in order that C(t) be quasimonotone nondecreasing with respect to K(B '), by 
Theorem C. 1 it is necessary and sufficient that BC(t)B _1 have its off-diagonal elements 
nonnegative. But BC(t)B _1 = A(t), which, by assumption, is not quasimonotone 
nondecreasing with respect toR”. □ 
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Hence, the problem is to find an appropriate cone with respect to which the 
function A(t) is quasimonotone nondecreasing, and for which the differential inequality 
D. 1 and all other required properties are preserved. The following case, where the 
inequality D. 1 is really an equality, is frequently found in applications. 

If the orbital derivative of a Vector Lyapunov function V € R “ satisfies V = AV 

where A is a stability matrix, then the following results give conditions under which the 
hypotheses of Theorem B.6 are satisfied, and the solution x(t) = 0 to equation A. 1 is a 
stable equilibrium. 

• • 

First, it is evident that if V = AV then B V < BAV for all nonsingular 

transformations B and all cones K. We also need the following result. 

Lemma D.3. If V(t, x) e R “ is locally Lipschitz in x with respect to R “ , then 

for any nonsingular matrix B e R nxn , Y = BV is locally Lipschitz in x with respect to the 
cone K(B). 

Proof. Evidently, V(t, x) e R” implies Y = BV e K(B). Since V is locally 
Lipschitz with respect to R “ , this implies there exists an LeR" such that 
L||x - y|| - (V(t, x) — V(t, y)) e R “ . Letting L = BL e K(B) yields 

L ||x - y|| - (Y(t, x) - Y(t, y» 6 K(B). □ 

We also require the result that the Lipschitz property is preserved through cone 
containment. 

Lemma D.4. If V(t, x) e P is locally Lipschitz in x with respect to the cone 
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P <z R", and P c= Q where QcR" is a cone, then V(t, x) is locally Lipschitz in x with 
respect to Q. 

Proof. If V(t, x) is locally Lipschitz in x with respect to P implies there exists 
an L 6 P such that L||x — y|| — (V(t, x) - V(t, y)) e P, then L e Q and 
L||x-y||-(V(t,x)-V(t,y))€Q. □ 

The following result appears as a remark in Lakshmikantham and Leela (1977). 
Lemma D.5. If P and Q are cones in R n with P <z Q, then x < y implies 

p 

x < y. 

Q 



Proof, x < y => y-xeP => y-xeQ =>x<y. □ 

P Q 



We now demonstrate the well-known result that if P and Q are cones in R" with 
PcQ, that the quasimonotonicity of a mapping A e R nxn with respect to Q does not 
imply the quasimonotonicity of A with respect to P. Let P = R* and let Q be the 



cone generated by the columns of B = 



1 0 
-3 1 



, hence PcQ. Let A = 



0 1 
-4 -6 



which is not quasimonotone nondecreasing with respect to P. Since B 1 AB = 



-3 1 

5 -3 



then by Theorem C. 1 A is quasimonotone nondecreasing with respect to Q. 

Similarly, if P cz Q are cones in R", the quasimonotonicity of a mapping with 
respect to P does not imply the quasimonotonicity of the mapping with respect to Q. 



Letting A = 



-2 2 
1 3 



and choosing P and Q as above, A is quasimonotone 
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nondecreasing with respect to P, but B 'AB = 



-8 

-32 



2 

9 



, so by Theorem C.l, A is not 



quasimonotone nondecreasing with respect to Q. 

Following Heikkila (1980) and Koksal and Fausett (1995) we present the 
following result. 

Proposition D . 6 . Suppose a Vector Lyapunov function V(t, x) e R ” satisfies 



V = A(t)V where A(t) is a stability matrix which is not quasimonotone nondecreasing 
with respect to R”. If the generalized left eigenvectors of A(t) are all contained in R " , 
there exists a transformation Y = B” 1 V and a cone K = R “ such that the hypotheses of 
Theorem B.6 are satisfied. 

Proof. The matrix A has as its Jordan canonical form J = B _1 AB, where the rows 
of B _1 are the generalized left eigenvectors of A and the columns of B are the 
generalized eigenvectors of A (see Horn and Johnson, 1991). By assumption, B - 1 is a 
nonnegative matrix, so K(B _1 ) cR”. Hence, the transformation Y = B -1 V yields the 

differential equation Y = JY. Then Y e K(B -1 ) is locally Lipschitz in x with respect to 
K(B -1 ) by Lemma D. 4, and with respect to R“ by Lemma D.5. Since J is 

quasimonotone nondecreasing with respect to R ” , and Y JY, we use Y as the 

vector Lyapunov function with K = R “ , satisfying the conditions of Theorem B.7. The 
only remaining detail is the measure used to bound Y, so the following remark completes 
the proof. □ 
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Remark D.7. Since PcQ implies Q' c; P*, we cannot, in general, find a 
| eQj =R“\{0} such that <4>,B -1 V) preserves the original measure 

n 

<4>o, V) = ^ Vi with <()o = [1, . - ., 1] T . This is because (j> = B<t> 0 is not generally in 
1=1 

Q * = R “ \{0}, since the nonnegativity of B' 1 does not imply the nonnegativity of B. 
However, had we used the equivalent measure max |V|| = max V;, then using as the new 

measure max |Yj|, since each Yj is a positive linear combination of the V;, preserves the 

i 

result of Theorem B. 6 by scaling a and b appropriately in equation B. 4. □ 

While we used the Jordan form to achieve quasimonotonicity with respect to R ° , 
this result may be extended to any case where S = B -1 AB, S has nonnegative off-diagonal 
elements, and the columns of B -1 are contained in R°. This follows from Theorem C.5 
of Heikkila (1980), who further shows that this construction is always possible in the case 
n = 2 if A has one eigenvector in the first quadrant and at most one of a i2 and a 2 i is 
negative (in this case it is necessary and sufficient that both eigenvectors be nonnegative). 
This is an application of Corollary C.5 to square systems with differential equalities. 

In the case where we do not begin with an equality in equation D. 1, this 

construction is generally not possible. Instead, given V < A(t)V, we now must find a 

R“ 

cone with respect to which A is quasimonotone nondecreasing, on which V is defined 
and locally Lipschitz in x, and on which the inequality D. 1 still holds. The result of 
Koksal and Lakshmikantham (1996) is a special case of this. 
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We define the following four properties for a linear comparison system for a vector 
Lyapunov function V and a cone K. 

(PI) V is defined on K; 

(P2) V is locally Lipschitz with respect to K; 

(P3) V | A(t)V; and 

(P4) A(t) is quasimonotone nondecreasing with respect to K. 

If we find a cone K = P for which properties (PI) through (P3) hold, then in 
order to find a cone K = Q for which properties (PI) through (P4) hold it is necessary 
that PcQ. This is because if a cone Q containing P with respect to which A(t) is 
quasimonotone nondecreasing exists, then by Lemmas D.4 and D.5, properties (PI) 
through (P3) hold for Q, and we can apply Theorem B.6. 

This is essentially what was done in Proposition D.2 where the vector Lyapunov 

function Y e K(B) cR' satisfied Y ^ ^ JY, and J was quasimonotone 

nondecreasing with respect to R“ . Letting K(B) = P and R" =Q satisfies properties 
(PI) through (P4). 

We first state an obvious result for finding such a cone, which may be useful in 
applications. The proof follows from the above remarks that PcQ. 

Proposition D.8. If V e P satisfies properties (PI) through (P3) for K = P, and 
if the cone K(B) generated by the generalized eigenvectors of A contains P, then 
properties (PI) through (P4) hold for K(B). □ 
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This generalizes to any transformation S = B -1 AB where S has nonnegative off- 
diagonal elements and K(B) contains P. Following Heikkila (1980) we next show that 
to find a cone containing R “ with respect to which A is quasimonotone nondecreasing, 
it is necessary and sufficient to find a cone contained in R° with respect to which A T is 
quasimonotone nondecreasing. We use Lemma n.A.3 which states that K(B) cR" if 
and only if K(B -1 ) dR”. Since it is presumably easier to find a cone contained in R ” 

(the matrix of its generators is nonnegative) than one containing R " (as we require) t hi s 
proposition is useful in applications. 

Proposition D.9. The matrix A is quasimonotone nondecreasing with respect to 
a cone K(B'')dR" if and only if A T is quasimonotone nondecreasing with respect to a 

cone K(B t ) cR". 

Proof. Let A T be quasimonotone nondecreasing with respect to a cone 
K(B t ) c R° . Then by Theorem C. 1, C = (B T ) -1 A T B T has all of its off-diagonal elements 

nonnegative, as does C 1 = BAB -1 . Since K(B T ) cR" implies K(B) cR”, then 
K(B -1 ) jR” and A is quasimonotone nondecreasing with respect to K(B ‘). 

Reversing the argument completes the proof. □ 

Hence we can apply Corollary C.4 and Theorem C.5 of Heikkila (1980) to the 
search for an appropriate cone, as the generalized eigenvectors of A T , or the generalized 
left eigenvectors of A, being contained in R“ is a sufficient condition for the existence 
of a cone for which property (P4) holds. 
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Because it is the case most frequently encountered in applications, we have 
assumed P = R" to establish properties (PI) through (P4). Evidently, if PcR" the 
above results still hold; however, if P 3 R ” then the construction is not as simple. 

Since for two cones P and Q, PcQ implies Q* c; P*, then if P R ” , instead 

»'▼’* -I (y *-p 

of seeking a set of vectors bj e R " such that (B ) A B has nonnegative off-diagonal 

elements, where the bi are the columns of B T , we must find vectors b, e P* to construct 
an appropriate cone. We therefore conclude the discussion of square linear comparison 
systems with the following extension of Proposition D.9. 

Theorem D.10. Given a vector Lyapunov function V and a cone PcR m for 
which properties (PI) through (P3) hold, in order to find a cone QdP with respect to 
which the matrix A is quasimonotone nondecreasing it is necessary and sufficient to find 
m independent vectors b, € P* such that if B T = [bi . . . b m ], then C = (B T ) _1 A T B T has 
nonnegative off-diagonal elements. 

Proof. If bi e P are such that C = (B )” A B has nonnegative off-diagonal 
elements, then so does C T = BAB -1 , and A is quasimonotone nondecreasing with respect 
to K(B _1 ) by Theorem C.l. Since K(B _l ) = K(B T )*, and since K(B T )eP* implies 
K(B t )* 3 P, then K(B _1 ) 3 P. Reversing the argument completes the proof. □ 

Proposition D.9, where P = R ” is a special case of the above result, since 
R “ = R “ *. Finding an invertible matrix whose coliunns are contained in P* is perhaps 
less difficult than finding one whose columns generate a cone containing P. The next 
result follows from Corollaries C.3 and C.4 of Heikkila (1980). 
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Corollary D.ll. For property (P4) to hold relative to some cone QdP, where P, 
Q, V, and A(t) are as above, it is necessary that at least one eigenvector of A T (t) belong 
to P and sufficient that A T (t) has n eigenvectors in P*, which do not depend on t. □ 
To show quasimonotonicity with respect to an appropriate cone, we only show 
such a cone exists. In an equivalent approach we could use the cone to transform the 
vector Lyapunov function and the comparison system as we did in Proposition D.2. This 
technique now works, since the cone used to make the transformation is contained in our 
original cone, and the inequality is preserved. The following corollary is an extension of 
Theorem 2.8.3 of Lakshmikantham, Matrosov, and Sivasundaram (1991). 

Corollary D.12. If a vector Lyapunov function V exists such that properties 
(PI) through (P3) hold for some cone P, and if a matrix B exists as in Theorem D. 10, 
then the vector Lyapunov function Y = (B T ) _1 V gives the same conclusion about the 
stability of the trivial solution to Equation A. 1. □ 

This approach has the advantage of directly yielding the stable manifold of the 
origin in the case of conditional stability, but in the usual case it involves an additional 
computational step. Since the ideas are equivalent, we continue the approach of showing 
such a cone exists. 

We continue our analysis of linear comparison systems considering the case where 
the system is rectangular instead of square. Let a nonnegative vector Lyapunov function 
V(t, x) be such that 

V: R+ x R m — » R ” , (D.5) 

• 

whose orbital derivative satisfies inequality D.4, V < A(t)V. Let A(t) e R nxn with 
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A = 



A A 

6 ••• o 



(D.6) 



where A e R pxp , A e R p * q , 6 is a qxp zero matrix, and 0 is a q x q zero matrix, 
with p + q = n. 



Such a comparison system is possible when the orbital derivative satisfies V j < 0 

for p < j < n. If we have equality in any of these components, or if V s = 0 for some j, 
this is a conserved quantity for our original system A. 1, and the theory of conservative 
systems to applies to our analysis. While it is necessary that a conserved quantity be non- 
constant on open sets, a component of a non-trivial vector Lyapunov function certainly 
satisfies this requirement, so this more powerful theory applies to the problem. Hence, we 
assume the orbital derivative satisfies inequalities in its last q components. 

This rectangular p x n system can be treated as a square system D.4 where A is 
given by equation D.6. Since we address such systems in Chapter IV, where we give a 
complete theory for reducible, square matrices, there is no need to consider the s mall er, 
rectangular system. 

We similarly dismiss the case where the matrix A in equation D.5 is of the form 



A = 



A 



A 



6 

0 
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where A e R qxp and the last q columns of A contain zeroes. A is not a stability 
matrix since cr(A) = o( A ) 'u {Ch, . . 0 q }, and A may not possess quasimonotonicity 

properties, even if A does. However, letting Yj = V; for 1 < i < p yields the square 

system Y < A Y, to which we may apply the theory for square comparison systems, 
regardless of the nature of A . Hence, we gain nothing by keeping the last q 
components of V, so we may discard them, leaving the imbedded p xp square system. 

We have now established the setting for the problem of finding a cone with respect 
to which a linear comparison system is quasimonotone nondecreasing. It is sufficient to 
consider square linear systems, and considering the transpose of the comparison system, it 
is sufficient to find a cone contained in the nonnegative orthant. The next two chapters 
present our solution to this problem 
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IV. THE QUASIMONOTONICITY OF A SQUARE, LINEAR OPERATOR 
WITH RESPECT TO A NONNEGATIVE CONE: THE REAL SPECTRUM 



For a matrix A e R nxn we address the problem “when is A quasimonotone 
nondecreasing with respect to a cone in R ” ?” (or equivalently, when is it essentially 
nonnegative under a nonnegative change of basis?). In Chapter HI we motivated this 
question with the application of determining stability in dynamical systems via the 
technique of cone- valued vector Lyapunov functions, and in Chapter VI we present 
further applications of this problem. 

Using Perron-Frobenius theory, Heikkila (1980) showed a necessary condition (the 
matrix A has a nonnegative first eigenvector) and a sufficient condition (the matrix A 
has all nonnegative eigenvectors) for such a cone to exist (see Corollary IH.C.2) where the 
first eigenvector is the nonnegative eigenvector xi associated with the (real) eigenvalue 
Xi of greatest real part. 

The next two chapters present our solution to this problem. Using constructive 
techniques, we bring together the necessary and sufficient conditions Heikkila gave for 
n > 2, and we address the complex spectrum for the first time in this setting. This chapter 
addresses the case where the spectrum of A is real, or g(A) c R, and Chapter V 
addresses the general spectrum. 

A. MATRICES WITH A POSITIVE FIRST EIGENVECTOR 

We begin with the case where the first eigenvector is positive, or xi > 0. We use a 
sequence of changes of basis Am+i = B A m B m with K(B m ) cR", which ensures the 
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change of basis is always nonnegative. The following le mma is the key to the 
construction. 



Lemma A.l. Let B = [x, l> 2 , . . ., b„] e R nxn be nonsingular with x e (R“ )° and 



B' 1 = 



4>i 

*1 






Then there exists a imit basis vector eteR" such that the matrix 



D 1 = 



(±)4>I 

(±)4>I 



has a nonnegative inverse D e R“ xn , where ( ± ) indicates an appropriate choice of sign 
for each <|)j. 

Proof. Since x e (R ” )° and <4>i, x) = 0 for i = 2 , . . n, then ±4); g R ” . 
Hence, each <f>j has components of both signs. Since each hyperplane <f> t contains 
x e (R ° )°, then each of these hyperplanes intersects at least (n - 1) of the coordinate 
hyperplanes e| in R“. (This is evident from the mixed signs in <{>; for i = 2, ..., n.) 

Therefore, there are (n - 1) hyperplanes which each intersect at least (n - 1) of the 
(ej) + = e r> R“ . By the pigeonhole principle, at least one (ej)^; intersects each <t> f . 

Let this unit basis vector be ek. Since we require dj) = 1 where d; is a column of D 
for i = 2, . . ., n, to ensure dj e (e k )^ we select the appropriate sign for each 4>i- 
Therefore, D _1 has as its inverse D = [di, d 2 , . . ., d n ], which is nonnegative. □ 
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We note that the kth row of D is ej , and that di = x/x*. 

We now construct a cone in R “ with respect to which a matrix A is 
quasimonotone non decreasing. We begin with the case where A is diagonalizable, and 
the first eigenvalue X.i is simple. In this case we can construct the matrix 



B' 1 = 









(A. la) 



where the <J> J are left eigenvectors of A and K((B -1 ) T ) dR" by Lemma A. 1. Then 
B = [xi, b 2 , . . ., b n ] is nonnegative, and 



A, = B”‘AB = 



X,, a 12 ••• a ln 
0 



(A. 2) 



If any a^ <0 for j = 2, . . ., n, then we use another change of basis of the form 



B,= 



1 ^ 



0 



C 

o 



(A. lb) 



where the are all nonnegative, and 



a 2 =b- 1 a 1 b i - 



A-i a 12 



“in 

0 



(A- 3) 
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Since aij = a ij + E,j(A,i — Aj), and Ai > A,j, we can make ay nonnegative, so A 2 is 
quasimonotone nondecreasing. 

In the case where A,i is not a simple eigenvalue, this construction still works, 
because for Aj = A.i, a ij = 0 in equation A.2. To show this, since aij = e^Abj and bj is 
in the eigenspace of Ay = Xi with [bj] k = 0 (since bj is orthogonal to all of the left 
eigenvectors c|)i where X, ^ A,), and to e^ ).' Hence a ij = Aibj = 0. 

When A is not diagonalizable, Koksal and Fausett (1996) extended Heikkila’s 
results using generalized eigenvectors. Following this, we use the “almost diagonal” form 
of A, or the canonical form of A with arbitrarily small 8 > 0 on the super- diagonal (see, 
for example, Horn and Johnson, 1991). Then, for negative a ij, aij = a ij + £j(A.i - Ay) or 
aij = a ij + £j(^i - Ay) - e£j_i, so that ay > 0 for arbitrarily small s. Since the only other 
nonzero off-diagonals of A 2 are 8 > 0, A 2 is again quasimonotone. 

The results of the previous two paragraphs combine in the case where Xi is the 
eigenvalue which makes A defective. We s ummari ze our results as follows. 

Theorem A.2. Let A e R nxn have a real spectrum and let a >j < 0 for some i^j. 
In order for A to be quasimonotone nondecreasing with respect to a nonnegative cone, it 
is necessary that the first eigenvector Xi > 0, and sufficient that xj > 0. □ 

The disadvantage of the above construction is that we must find all of the 
generalized left eigenvectors of the matrix A in to construct the matrix A 2 (and the cone 
K(BB t )). However, since we only need to know whether such a cone exists for stability 
results, it may be sufficient to compute only the first eigenvector. 
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Since in the case where xi has some negative component we know no cone exists, 
and when X] is positive a cone always exists, the only remaining case for a real spectrum 
is when xi e <9R ” . 

B. REDUCIBLE MATRICES WITH A NONNEGATIVE FIRST 
EIGENVECTOR 



This section addresses the case where xi e 5R " , or when xi > 0 has one or 
more components equal to zero. We assume xi e (RJ)°© 0 q , where 0 q indicates the 
zero vector in R q , and p + q = n. Because Xi e 3R ” we cannot necessarily apply 
Lemma A. 1 to the left eigenvectors of A to obtain a nonnegative change of basis. 

Since under a (suitably ordered) nonnegative change of basis B, with xi e K(B), 
Xi i-» B _1 xi e R \ © 0 q . From this and Theorem n.B.2, it follows that if Ai = B -1 AB = 



L n 



12 



l 22 



is quasimonotone, then if A n is irreducible, it is no larger than p x p. 



Furthermore, if A 22 is irreducible, then it can be put into the upper triangular form of 



equation A. 3 via a transformation B = 



I 0 



0 B 



-i 
22 J 



, where B 22 has the form of equation 



A. 1. Since this transformation does not effect the nonnegativity of A 12 , in this case the 
quasimonotone matrix Ai = B _1 AB has a reduced form with a q x p zero block. From 
this follows our first result about reducible matrices, but we first clarify the particular form 
the reducible matrix must have. 



For a reducible matrix A with a first eigenvector x t e 3R” , when A is 



permuted so that xi has all ofits zeros in its last q components and A has a q x p 
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lower-left zero block, we define A as reduced. We treat matrices which are not reduced 
when xi is ordered as above with the irreducible case. 



Theorem B.l. A reduced matrix A = 



An ^12 



A 

A 



22 . 



is quasimonotone with 



respect to a nonnegative cone only if Ah and A22 are quasimonotone nondecreasing 
with respect to nonnegative cones. 

Proof. The previous remarks address the case where the quasimonotone form 



B'AB = A a = 



<< 0 
1 


A + 
^12 

A2 9 


is reduced. Consider the matrix A q = 


< ° 


a; 

A 2 q 2 


K 

A 23 




22 j 




0 


0 


A 33 



where A ? are quasimonotone and irreducible, and A* > 0 . Let the diagonal blocks 
have dimensions r x r, s x s, and t x t respectively, where r < P, t < q, r + s > p, and 
s + 1 > q. Since the largest eigenvalue Hi of A22 is either an eigenvalue of A\ 2 or 
A , then A<, has associated with this eigenvalue an eigenvector which is nonnegative in 
the last s + 1 components. Hence, A22 has a nonnegative eigenvector associated with 
Hi. The case where this eigenvector is positive is the case of A q above, and when A22 
is reduced, we can repeat this argument. Since the result of this theorem is trivial in the 
case p = 1 (which is the only case required to prove Theorem C.2 in the next section) 
then the case where A22 is not reduced follows from that theorem. □ 

We now construct a cone when A is reduced, with An e R pxp and A22 e R q *q 
Assume that a(An) = (>n > X2 > . . . > Xp} and a(A 2 2) = (hi > H2 ^ ^ Hq} where 

A-i > pi and A22 has a positive first eigenvector. 
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First let B 0 = 



B„ 0 
0 B 



22 J 



, where Bn and B 22 have the form of equation A. 1, so 



Aj = B 0 'AB 0 = 



X,j a 2 -- 
0 

0 



a. 



X 



1 12 



M-l P 2 *" Pq 
0 



has upper-left and lower-right blocks with the form of equation A.3. While this 
transformation does not leave A 12 unchanged, we can now let A 12 be arbitrary 



Next let Bi = 



I B 



12 



0 I 



where Bi 2 = 



with y; > 0 and 



A 2 = B, A,Bj . This only changes the first row of A 12 , where 
a i,p+i i-» a i jP+ i + yi(Xi - Hi) = ai, p+ i, and a ^ i-» a i >p+ j - yiPj + y^i - Hj) = ay for 



j = 2, . . q. Since Xi > p.j, ay can be made arbitrarily large. 



Now let B 2 = 



Dn 0 
0 I 



, where Dn = 



1 

-5, 



-5 p 0 



0 



with 8| > 0, and 



A 3 = Yi 2 A 2 B 2 . While B 2 is not nonnegative, we show that the product BoBiB 2 still is, 
but first we examine the effect of this transformation. The off-diagonals in the first row of 
A 2 are unchanged under B 2 , as are the last q rows. For i = 2, . . ., p, j = 1, . . ., q, 
a,,p+j 1 — > aj,p+j + 5, ay, and since ay is arbitrarily large, aj <p +j can be made nonnegative for 
arbitrarily small 8;. Furthermore, for i = 2, . . ., p; j = 2, . . ., p; i * j, ay i-» SjCtj > 0. 
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Finally, for i = 2, . . p; an 8j[Li - A,j - ^]a k 5 k ], and since X\ > ^ and the 8 k can 



k=2 



be made arbitrarily small, then an > 0 as well. Hence, A 3 is quasimonotone. 



We now show that B = B 0 BiB 2 is nonnegative. Since B = 



B U D 11 B 11 B 12 



0 B 



22 



we only consider BuDn. If we permute An so that when constructing Bn, ek = ei in 



Lemma A. 1, then Bn = 



1 0 

: 2 b 22 "* 7p 



X P b p2 



, where b;j > 0 and x; > 0, since An has a 



positive first eigenvector. Only the first column of Bn changes under BnDn, and 



Xj h-9> xj - 7,8, b„ >0 since x, > 0 and the 8j are arbitrarily small. Therefore, in this 

j=2 



case, a nonnegative cone always exists. 

We now show we can relax the assumption that An and A 2 2 have first 
eigenvectors which are strictly positive. If An and A 22 are each quasimonotone 
nondecreasing with respect to nonnegative cones, then by the above construction, each 
can be put in a quasimonotone form with arbitrarily small nonzero elements below the 
main diagonal. In this case, the above construction still works for the matrix A. We 
summarize our results. 



Theorem B.2. For a reduced matrix A = 



An ^12 



0 A 



22 , 



with real spectrum and 



first eigenvector Xi e (R p + )° © 0 q to be quasimonotone nondecreasing with respect to a 



nonnegative cone, it is necessary that An and A 22 are quasimonotone nondecreasing 
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with respect to nonnegative cones, and sufficient that Xi > pi, where Xi and pi are the 
first eigenvalues of A n and A22 respectively. □ 

The only reduced matrices which w'e have not included in this theorem are those 
where pi = X\. If both An and A22 are in upper-triangular quasimonotone form with 
Pi = A,i, this construction may or may not be possible, depending on the sub-matrix Ai 2 . 



The following example with n = 2 illustrates this. When X > p, A = 






M2 



is quasimonotone with respect to K(B) where B = 



1 

0 



y 

1 



for sufficiently large y, 



regardless of the sign of a^. Also, if an > 0 , then A is quasimonotone with respect to 
K(I), regardless of X and p. However, if p > X and a 12 < 0 , then A is not 
quasimonotone with respect to any nonnegative cone. (The case p > X violates the 
necessary condition of the first eigenvector being nonnegative.) Furthermore, if an > 0 , A 
has two nonnegative generalized eigenvectors, but if an < 0 , A has a generalized 
eigenvector x 2 ?R 2 + . 

Using the same construction as we did when A-i > pi, we formulate the following 
sufficient condition for a cone to exist, which we do not describe in terms of generalized 
eigenvectors, but instead in terms of the sub-matrix An- 



Corollary B. 3 . For a reduced matrix A = 



A 

z M2 

0 in . 22 



A 

A 



with real spectrum and 



first eigenvector Xi e (R p + )° © 0 q to be quasimonotone nondecreasing with respect to a 
nonnegative cone, it is necessary that An and A22 are each quasimonotone with respect 
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to nonnegative cones. Furthermore, if An and A22 are in upper-triangular 
quasimonotone form with an = ^i and ap+i 7 p+i = pi = A,i, where pi is a simple 
eigenvalue of A 22 , then it is sufficient that ai,p+i > 0 for i = 1 , . . ., p. □ 

We can formulate similar sufficient conditions from this construction if pi is not a 
simple eigenvalue of A 22 . Given the non-uniqueness of generalized eigenvectors, we have 
not produced a condition in this case that is based on generalized eigenvectors, although in 
the case where n = 2 the sufficient condition in Corollary B.3 is also necessary, and the 
requirement for both generalized eigenvectors to be in R + is evident. 



C. IRREDUCIBLE MATRICES WITH A NONNEGATIVE FIRST 
EIGENVECTOR 



This section addresses the case where xi e dR “ , but A is either irreducible or 



not reduced. Consider the matrices A = 


'7 

4 


-2 

1 


— 7" 
-9 


and A = 


"7 

6 


-2 

-1 


-7' 

-9 




-1 


.1 


1 




-1 


.1 


3 



, which 



have = 5, « 2.894, X 3 » 1.106, and xi = [1, 1, 0] T A is quasimonotone 



nondecreasing with respect to K(B) where B = 



1 2 15 

1 3 17 

0 .05 1 



, yet A is not 



quasimonotone nondecreasing with respect to any nonnegative cone. 

From the discussion preceeding Theorem B. 1, if A has a quasimonotone form via 
a nonnegative change of basis, then some such matrix A, is reduced with a lower-left zero 



block. Furthermore, if A q 





and An has a positive first eigenvector (for 
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example, if it is irreducible) then we can put An in the form of equation A.3 using a 
transformation of the form of equation A. 1. If Xi is simple, we can use the techniques of 
Section B for reduced matrices to ensure Ai 2 is nonnegative. This proves the following. 

Lemma C.l. If a matrix A with real spectrum and first eigenvector 
Xi € (RJ)° © 0 q is quasimonotone nondecreasing with respect to a nonnegative cone, 
then it is quasimonotone nondecreasing with respect to a nonnegative cone for which xi 
is an extremal. □ 

We cannot generally construct a nonnegative cone with respect to which A is 
quasimonotone as we did in Sections A and B (since the example in this section shows 
that such a cone does not always exist) but we can reduce the dimension of the problem by 



f(l) O' 

one by deflating the matrix with the first eigenvector as follows. Let Bo = 

_ x i l . 

where I is the identity matrix of dimension n - 1, and Xi = [1, £ 2 , • • £ P , 0, . . .,0] T , so 



Aj = B-'AB 0 = 



X 


p 




1 


a lj~ 


1 

o 


V 




o 
1 


1 

<N 

M 

< 



where a^ = ajj - aij£;, and for i > p, a ;j = aij. 



Since Theorem B. 1 is trivial in the case p = 1, Ai is quasimonotone with respect 
to a nonnegative cone only if A 22 is, so this is a necessary condition for A to be 
quasimonotone with respect to a nonnegative cone. Furthermore, since Xi > X; for i = 2, 
. . ., n (i.e., X; is an eigenvalue of A 22 ) then using a transformation Bi of the form of 
equation A. lb, we can make A 2 = B^AjBj quasimonotone as well. 

The cone K(B 0 ) is a subcone of R p + © R q , where we replaced one of the first p 
basis vectors with xi. The transformation we showed replaced ei, but the choice of 
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which basis vector to replace can be varied by permuting the first p components of xj. It 
is only necessary that one such transformation produce a matrix A22 which is 
quasimonotone with respect to a nonnegative cone in order for the above construction to 
produce a cone for the matrix A. We summarize these results as follows. 

Theorem C. 2 . Let A e R n<n have a real spectrum with a first eigenvector 
Xi e (RP)°© 0 q , and let A be irreducible (or not reduced). Then for A to be 
quasimonotone nondecreasing with respect to a nonnegative cone it is necessary that for 



B 0 



( 1 ) 0 

.*1 I 



for some permutation Xj e (R^ )° © 0 q of xi, and for 



A, - B 0 'AB 0 - 



A,j a j j 

0 A 



22 j 



, that the matrix A22 be quasimonotone nondecreasing with 



respect to a nonnegative cone. In this case, it is sufficient that A, t be a simple eigenvalue 
of A. □ 



When k\ is not simple the problem is the same as the one discussed in the reduced 



case, where, for example, A = 



X a 



12 



0 X 



, and the quasimonotonicity depends on the sign 



of ai 2 . 

The choice of which permutation of Xi to select is not obvious, but there are only 
p such choices that lead to distinct cones. One solution is to compute the second 
eigenvector x 2 of A, associated with X2, which is the eigenpair required for A22 to be 
quasimonotone. Since the last n - 1 components of B~'x2 are the first eigenvector of 
A22, then finding a B " 1 such that B" l x 2 > 0 is necessary, and that the last n - 1 
components of B _ 1 x 2 > 0 is sufficient, to determine that such a cone exists. 
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V. THE QUASIMONOTONICITY OF A SQUARE, LINEAR OPERATOR WITH 
RESPECT TO A NONNEGATIVE CONE: THE GENERAL SPECTRUM 



A. IRREDUCIBLE, ESSENTIALLY NONNEGATIVE MATRICES. 

In this chapter we discuss the problem of finding a cone with respect to which the 
matrix A is quasimonotone nondecreasing in the case where A has a general (complex) 
spectrum, or when g(A) is not strictly real. Because A is real, its complex eigenvalues 
and eigenvectors occur in conjugate pairs; however, unlike with the real spectrum, the 
real-diagonal form of A is not quasimonotone (it has diagonal blocks which necessarily 
have a negative off-diagonal element) so we can no longer use the diagonal form. 

Additionally, a matrix A with a general spectrum is not always similar to a 
quasimonotone matrix, much less through a nonnegative change of basis. The problem of 
determining when a given set of complex numbers is the spectrum of a nonnegative matrix 
is known as the nonnegative inverse eigenvalue problem, and we discuss it in the next 
section. Since a matrix with a real spectrum is always similar to a quasimonotone matrix 
(its diagonal, or Jordan canonical form) this was not an issue in Chapter IV, and the only 
task was to produce a nonnegative change of basis. However, in this chapter, the 
nonnegative inverse eigenvalue problem plays a significant role in our solution. 

We have some a m ount of flexibility over the nonnegative inverse eigenvalue 
problem in that the matrices we seek are essentially nonnegative, and we can make them 
nonnegative by shifting the diagonal (and the spectrum) by a positive quantity r. Since 
B _1 AB is essentially nonnegative if and only if B 'AB + rl is nonnegative for sufficiently 
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large r, and since B _I (A + rI)B = B _1 AB + rl, for this problem to have a solution the 
nonnegative inverse eigenvalue problem must have a solution for the set o(A) + r. 
We now show that the matrices we need to consider have an irreducible, 

essentially nonnegative form. Otherwise, if B -1 AB = A, = 



An A j2 



L 22 



is essentially 



nonnegative, then the shifted spectra a(An) + r and o(A 22 ) + r must both solve the 
nonnegative inverse eigenvalue problem. In this case we transform A ]2 using the 
construction in Chapter IV, although if An and A 22 are not reducible, we cannot take 
advantage of the upper-triangular structure they had in the case of the real spectrum We 
use this, and the methods of Chapter IV, to discard a portion of the real part of o(A), and 
we consider only that part which is required to form a solution to the inverse eigenvalue 



problem For example, if a(A 22 ) c R, by using transformations B = 



I 0 
0 B 



22 . 



where 



B 22 is formed using the techniques of Chapter IV, we only need the theory for the 
complex spectrum to treat An. 

Since the final quasimonotone form is irreducible, it has a positive first 
eigenvector, so A having a positive first eigenvector is a necessary condition for this 
problem to have a solution. (Or as in A<, above, Xi must have at least enough positive 
components to correspond to the dimension of An.) Furthermore, from Frobenius 
(1912) the first eigenvalue must be simple or the quasimonotone form is reducible. 

Since the (shifted) spectrum of A must solve the nonnegative inverse eigenvalue 
problem for A to be quasimonotone nondecreasing with respect to a nonnegative cone, 
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and it is necessary that the first eigenvector of A be positive when the quasimonotone 
form, Aq, is irreducible. We now show that these conditions are sufficient. 

Without loss of generality, assume A, is not cyclic (otherwise we can further shift 
its spectrum so that it is non-cyclic). Let D be the real-diagonal canonical form of both 
A and Aq, so D = B AB 0 = B AqBi. Not only are the spectra of A and Aq the same, 
but the matrices are similar, so that the eigenvalues have the same geometric multiplicity 
as well. Hence, when we require the shifted spectrum of A to solve the nonnegative 
inverse eigenvalue problem, this includes the geometric multiplicity of the eigenvalues. 
Here, Bo = [xi, . . ., xj and Bi = [b b . . ., b n ] where the Xj and b; are the real parts of 
the (generalized) eigenvectors of A and Aq respectively, and Xi and bi are the 
positive first eigenvectors associated with the simple eigenvalue A-i of greatest real part. 

While Aq = BiB g 1 ABoB , it is not generally true that BoB > 0, as required. 
However, we further change bases using A” as n — » oo without changing Aq to make 
the change of basis nonnegative. Let B = BoB j -1 A° , so 

Aq = B _1 AB = (A " ) _1 B iB AB 0 B A ” . Since Aq is irreducible, non-cyclic, and has a 
simple first eigenvalue, A” — > [bi, . . ., bi] as n -> co. Therefore, as ,n — > oo, 

= [xi, . . ., xi] > 0, so A is quasimonotone 

nondecreasing with respect to a nonnegative cone. 

Hence, given a positive first eigenvector, we have reduced the problem to the 
nonnegative inverse eigenvalue problem,. We summarize our result as follows. 



B — » BoB / [bi, . . ., bi] - Bo 



67 



Theorem A.l. For a matrix A e R nxn with complex eigenvalues to be similar to 
an irreducible, essentially nonnegative matrix by a nonnegative change of basis, it is 
necessary and sufficient that A have a positive first eigenvector (associated with the 
simple real eigenvalue of greatest real part) and that the shifted spectrum of A solve the 
nonnegative inverse eigenvalue problem for an irreducible matrix. □ 

In the case where the final quasimonotone form is reducible, A 2 2 as well as An 
may have a non-real spectrum, in which case we can apply the above theorem to each of 
them individually, and then use the techniques of Chapter IV to make Ai 2 nonnegative. In 
the case where a(A 22 ) is strictly real we can use the techniques of Chapter IV to “throw 
out” this real part of the spectrum, and construct An using the above techniques. The 
following example shows, however, that we cannot simply “throw out” all but Xi from 
the real part of the spectrum. 



The matrix 



'0 1 
0 0 
0 0 
1 0 



0 

1 

0 

0 



0 

0 

1 

0 



is a 4x4 circulant matrix (each row is a shift of 



the previous one; see Davis, 1979), which is a type of cyclic matrix (its eigenvalues are the 
four roots of unity). Since a(A) = (±1, ±i}, this set is a solution to the nonnegative 
inverse eigenvalue problem. However, we cannot discard the eigenvalue A, 4 = -1, since 
the set c = { 1, i, -i} is not the spectrum of any nonnegative (or quasimonotone) matrix. 
Furthermore, the left eigenvector associated with A , 4 has two components of each sign, so 
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that the transformation required to throw out B 1 = 

nonnegative inverse B. 

When a matrix has a reducible quasimonotone form under a nonnegative change of 
basis we cannot discard all but the essential real part of the spectrum prior to constructing 

A<,. Using transformations of the form B _1 = 

left eigenvectors, requires B to be nonnegative, which is not true in general. However, in 
the cases where this construction can be done a priori, then Theorem A. 1 gives sufficient 
conditions for our problem to have a solution. 

To complete this topic, we present a brief discussion of the nonnegative inverse 
eigenvalue problem. 

B. THE NONNEGATIVE INVERSE EIGENVALUE PROBLEM 

There are many types of inverse eigenvalue problems. Given a set of numbers and 
a structure, they ask when that set of numbers is the spectrum of a matrix with the given 
structure. For a thorough survey of inverse eigenvalue problems, see Chu (1998). The 
structure we require is essential nonnegativity, or nonnegativity for the shifted spectrum. 
The nonnegative inverse eigenvalue problem is one of the classic unsolved problems from 
the theory of linear algebra (see, for example, Horn and Johnson, 1991) although some 



c 



for some collection of n - p 






0 



♦I 



cannot have a 
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sufficient and some necessary conditions have been given for a set of complex numbers to 
be the spectrum of a nonnegative matrix. 

The nonnegative inverse eigenvalue problem dates back to Perron and Frobenius, 
and has been approached in numerous forms and variants since then. Necessary conditions 
for real spectra were given by Suleimanova (1949). His results were improved by 
Karpelevich (1951), Perfect (1952, 1953, 1955), Suleimanova (1965), and Fiedler (1972). 
Fiedler (1974) extended these results with the additional requirement that the matrix be 
symmetric, and Friedland and Karlin (1975) solved a variant of the problem known as the 
nonnegative additive inverse eigenvalue problem. Most of the results for real spectra are 
refinements of Suleimanova’s necessary condition that the absolute value of the sum of the 
negative eigenvalues may not exceed Xi > 0. Results of this type are always obtainable 
for essentially nonnegative matrices since we are allowed to shift the spectrum, and as 
mentioned earlier, the diagonal matrix of any real spectrum is essentially nonnegative. 

For complex (not strictly real) spectra, the problem is more difficult. Because of 
the close relation to nonnegative matrices (see von Mises, 1931) it has been addressed in 
various forms for nonnegative stochastic matrices, starting with Frechet (1933), 
Romanovsky (1936), Dmitriev and D ynkin (1945), Taussky (1948), and Brauer (1952). 
Brauer presents a summary of the previous results, in addition to his own results and 
methods of proof. Friedland and Melkman (1979) resolve the question in the case of 
nonnegative Jacobi matrices, and numerous bounds for the spectra of nonnegative 
matrices exist (see, for example, Ostrowski and Schneider, 1960) as well as results on the 
first eigenvector of nonnegative matrices (see, for example. Mine, 1970; or Ashley, 1987), 
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but these latter results require a priori knowledge of the matrix. Additionally, Friedland 
(1977) extended his previous results for real spectra to complex spectra for the additive 
and multiplicative nonnegative inverse eigenvalue problem 

Of the more recent results on the nonnegative inverse eigenvalue problem for 
complex spectra, as it applies to our problem, we mention those which provide sufficient 
conditions that may be used by the analyst, for example, to draw a stability conclusion in 
the method of vector Lyapunov functions. Kellogg ( 1971) addressed the problem of when 
a set of complex numbers is the spectrum of an essentially positive matrix (hence, for our 
application, certainly an irreducible essentially nonnegative matrix). He found a sufficient 
condition to be Re (Lj) - Re(L;+i) > V3 Im(Li+i) when the eigenvalues are ordered by 
their real parts, where A-i is the first (real) eigenvalue, and the X, are the eigenvalues 
with positive imaginary parts for i > 1. The other real eigenvalues do not affect this 
condition. Certainly, the condition is not necessary, as the set { V2 + e, ±i} is the 
spectrum of a positive matrix for all 8 > 0. 

Friedland (1978) addressed the problem for eventually nonnegative matrices 

D 

(where A k is nonnegative for all k > M for some M >1) and showed that if ^ A. 1 ’ >0 

i=i 

for some p e { 1, 2, . . .} and if X\ is the only positive Lj for i = 1, . . ., n, then A is 
eventually nonnegative. This supports the conjecture that if k = 1, then {Li} is the 
spectrum of a nonnegative matrix, although no proof has been found yet. Furthermore, if 
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A > 0, A k > 0, it is necessary that ^ A. k >0. Similarly, Lowey and London (1978) gave 

i=i 

n n 

the further necessary condition that ( ^ A, k ) m < n m_1 ^ for all k, m e { 1, 2, . . . }. 

i=l i=l 

Boyle, et. al., (1991) characterized via symbohc dynamics the sets of complex 
numbers which are the nonzero portion of the spectrum of a nonnegative matrix; however, 
their result is conjecture. Their results are important in that they are presented for 
primitive matrices (those matrices A > 0 for which A k > 0 for some k). For the 
nonnegative inverse eigenvalue problem, this is sufficient as the generalization to 
irreducible or nonnegative matrices is easily determined. More significant, however, is 
that for our stability application, since a primitive matrix is a special case of an irreducible 
nonnegative matrix, this is all we need to apply the techniques of the previous section. 

Further results, for example the sufficient conditions of T. J. Laffey or those of 
Koltracht, Newman, and Xiao (1993) (where Boyle’s conjecture is shown to be true for 
n < 5), are discussed in Berman and Plemmons (1994). There are probably enough 
necessary and sufficient conditions in the literature that most spectra encountered can be 
determined to be solutions of the nonnegative inverse eigenvalue problem (or not), but as 
long as the problem is unsolved, the pathological case cannot be discounted. 
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VI. FURTHER APPLICATIONS 



This chapter presents additional applications that can benefit from the results of 
this dissertation. 

Boyd, et. al., (1994) present the problem of positive orthant stabilizability, which 
is a special case of the “hit and hold” problem from control theory (see Berman, Neumann, 
and Stem, 1991). Given a linear differential system x' = Ax where A e R nxn , when does 
x(0) eR“ imply x(t) e R“ for all t > 0, and x(t) — > 0 as t — » co? It is a necessary and 
sufficient condition that there exist a diagonal matrix P > 0 such that PA T + AP < 0, and 
that A be essentially nonnegative. 

Ohta, Maeda, and Kodama (1984) discuss a s imil ar problem where the state 
variable is not required to be nonnegative, in which case the problem of finding an 
essentially nonnegative matrix is the nonnegative inverse eigenvalue problem for the 
shifted spectrum They also present the positive realization problem, which is similar to 
the problem of positive orthant stabilizability, and for which essentially nonnegative 
matrices play an important role. 

In the problem of positive orthant stabilizability, we cannot arbitrarily shift the 
spectrum of a linear dynamic system without changing the problem However, if the 
matrix A is not essentially nonnegative, but its spectrum satisfies the necessary 
conditions, then our techniques can find a cone contained in the positive orthant which can 
be stabilized under the original system Our techniques can also be used to find a cone 
containing the positive orthant which can be stabilized under the original system 
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A related problem from control theory in Stem (1980) is that of asymptotic 
holdability, where a trajectory must be held arbitrarily close to a linear subset of R n (not 
necessarily the nonnegative orthant) for any finite time interval. This leads to positively 
invariant cones, which Stem (1982) characterized for cones other than R” . These are 
cones K such that when x(t) is the solution to x' = Ax, then x(0) e K implies x(t) 6 K 
for all t > 0. A sufficient condition in this case is the essential nonnegativity of A with 
respect to K. 

The hit and hold problem considers a positively invariant simplicial cone K in R", 
and determines the set of initial values X A (K) in R n which eventually reach K (and 
hence remain there due to the positive invariance of K) under x' = Ax. Neumann and 
Stem (1985) show that X A (K) is itself a cone, and they characterize it for diagonalizable 
systems A. Berman, Neumann, and Stem (1986) extend these results to 
nondiagonalizable systems A. These results require the spectrum of A to be real, and a 
result from Berman, Neumann, and Stem (1991) requires the matrix A to be essentially 
nonnegative as a condition of the theorem. 

In the case where A is not essentially nonnegative, the hit and hold problem has 
not been solved. However, we use our techniques to find a nonnegative cone Ki with 
respect to which A is essentially nonnegative, and the theorem of Berman et. al., can be 
applied to a smaller cone Ki. If X a (Ki) 3 K, then this is a solution to the hit and hold 
problem for K_ 

These results have applications to economics (see, for example, Sierksma, 1979), 
engineering, and biology, where frequently state variables are nonnegative and the 
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essential nonnegativity of the differential system is often required. If a problem is defined 
on the nonnegative orthant and the operator is not essentially nonnegative, our solutions 
may be used to find a sub-cone of R “ for which the problem has the required properties, 
and to which the results from control theory may be applied. 

Bitsoris (1988) studied the hit and hold problem for discrete-time dynamical 
systems x^+i = Ax^, and showed that in addition to invariant cones, invariant polyhedral 
sets also arise, and these results could be applied to the unsolved problem of determining 
the spectral properties of systems possessing positively invariant polyhedral cones. For 
the continuous-time case, Castelan and Hennet (1993) again related the existence of 
positively invariant polyhedral sets to the problem of finding an essentially nonnegative 
matrix. 

Since we have a technique for finding essentially nonnegative matrices with respect 
to smaller cones (or larger cones if we consider the transpose of the original system) then 
depending on the requirements of the application we also have a useful technique for 
finding positively invariant polyhedral sets in dynamical systems. 

The requirement that a change of basis be nonnegative is a natural one for many of 
these applications, although the techniques we developed can be extended to cones 
outside of the nonnegative orthant (where for the real spectrum, the problem of seeking an 
essentially nonnegative matrix becomes trivial, and for the general spectrum is simply the 
nonnegative inverse eigenvalue problem for the shifted spectrum). Furthermore, our 
results for complex spectra may be used to extend current results which exist only for 
matrices with real spectra. 
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Stem (1982) shows that when the general cone K = R ° , the problem of positively 
invariant cones has its solution in the theory of M-matrices. An M-matrix (after 
Minkowski) is a square matrix A = rl - B where B > 0 and r > p(B) (or r > p(B) for 
nonsingular M-matrices). This is a subset of the class of essentially nonpositive matrices 
with nonnegative diagonals. M-matrices occur frequently in problems from the biological, 
physical, and social sciences (see, for example, Plemmons, 1977) as well as from the 
mathematical sciences, and Varga (1976) points out that there is an infrequent exchange of 
results between these disciplines. 

Stem (1981) shows that while a matrix which is cone invariant under the 
differential system x' = Ax, x(0) e K is more general than the requirements for A to be 
an M-matrix. However, the main properties which hold for M-matrices also hold for A, 
which he calls a generalized M-matrix. Numerous classifications of M-matrices have been 
shown (see, for example, Rothblum, 1979; or Stem and Tsatsumeros, 1987) and most of 
these are related to the cone invariance of A or the exponential nonnegativity of -A. 

In terms of other applications of M-matrices, Araki (1975) shows a direct 
application to the method of vector Lyapunov functions (similar to what we presented in 
Chapter III), and Kielson and Styan (1973) relate them to the theory of Markov chains 
which they use to establish results concerning the nonnegative inverse eigenvalue problem 
for complex spectra. 

The negative of an M-matrix is an essentially nonnegative matrix, but the problem 
we present has a less obvious but more important connection to the theory of M-matrices, 
although our solution technique does not provide an immediate solution. The inverse M- 
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matrix problem (see Willoughby, 1977) asks when the inverse of a positive matrix is an M- 
matrix, and it is also mentioned that nonnegative matrices which are not strictly positive 
may also have M-matrices as their inverses, but the zeros may only be those associated 
with reducibility. Since the spectra of M-matrices have certain properties, not every 
positive matrix solves the inverse M-matrix problem, but as our problem allows a shift in 
the spectrum, an irreducible solution to our problem may be a way to obtain an M-matrix 
(or its inverse) which might have use in the numerous applications of M-matrices. 

A more general class of matrices than M-matrices is those with nonpositive off- 
diagonal entries (see Berman, Varga, and Ward, 1978), or essentially nonpositive 
matrices. This class has an important subclass where all of the principal min ors are 
positive (see Fiedler and Ptak, 1962) which is also a more general class than M-matrices. 
Again, the relationship of these classes of matrices to essentially nonnegative matrices is 
evident, particularly if we are allowed to shift the spectrum. While the techniques we 
present do not easily adjust to find a matrix with nonpositive off-diagonal elements, it may 
be possible to develop a similar technique for essentially nonpositive matrices. 

Further classes of matrices studied in relation to the problem of when a matrix is 
positive on a cone include the cross-positive matrices , which are those matrices A for 
which xeK and y e K* are such that (x, y) = 0, then (x. Ay) >0 for a proper, 
simplicial cone K. These matrices are discussed in Schneider and Vidyasagar (1970) and 
Tam (1975), and the former draws the relationship to essentially nonnegative matrices. 
This relationship is through the exponential of a matrix, where a matrix A is cross- 
positive on K if and only if e A is positive on K. This is true if and only if A is 
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essentially nonnegative on K. While the paper shows the case of K being the 
nonnegative orthant, the change of bases we constructed can make the result more 
general, either in terms of smaller or larger cones. The class of cross-positive matrices 
includes the copositive matrices, which are symmetric and for which K = R“ (see 
Haynsworth and Hoffinan, 1969). These matrices are also essentially nonnegative, and in 
addition to the applications previously mentioned, they are important in the field of 
mathematical programming. 

The question of when a linear differential operator is quasimonotone with respect 
to a nonnegative cone is related to the question of when a linear transformation has an 
invariant cone. Birkoff (1967) studied this problem and used the results to give a 
constructive proof of Perron’s theorems. Following Birkoff Vandergraft (1968) 
characterized the spectral properties of matrices with invariant cones. These results were 
extended by Barker and Turner (1973) and were related to the nonnegative inverse 
eigenvalue problem. These ideas have been used in algebraic Perron-Frobenius theory 
(see Barker and Schneider, 1975) and have been studied geometrically to extend the ideas 
to non-simplicial cones (see Barker, 1973). An algebraic extension ofLowey and 
Schneider (1975) further characterizes matrices with invariant cones. These techniques 
were applied by Vandergraft (1972) to establish the convergence of some splitting 
methods for solving systems of linear equations. For example, both the Jacobi and Gauss- 
Seidel methods converge for essentially nonpositive matrices with positive inverses (see 
Varga, 1962), although to apply our techniques to this problem may be more challenging 
numerically then solving the original problem Ax = b. 
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The question of when a matrix is quasimonotone nondecreasing (or nonin creasing) 
with respect to a cone is not new, although the requirement that the cone be nonnegative 
does not appeared frequently in recent applications. In applications where the state 
variables are required to be nonnegative (for example, in population dynamics), and when 
the original problem is defined on the nonnegative orthant, this requirement is important. 
These applications, in addition to the one of cone-valued vector Lyapunov functions, 
demonstrate the usefulness of the results of this dissertation. 
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VII. CONCLUSIONS AND FUTURE RESEARCH 



This dissertation asks, and answers, the question of when a matrix 
A e R™" is quasimonotone nondecreasing with respect to a nonnegative cone. Chapter 
IE motivates this question with the theory of cone-valued vector Lyapunov functions from 
stability theory in dynamical systems. Chapter VI presents other applications and related 
problems where the results can be used. We summarize these results as follows. 

A matrix A e R nxn is quasimonotone nondecreasing with respect to a nonnegative 
cone if the following conditions hold. 

(1) (Heikkila, 1980). It is necessary that the matrix A have a real first eigenvalue Xi 
such that X\ > Re(Aj) for i = 2, . . ., n. Furthermore, it is necessary that associated with 
X-i be a first eigenvector Xi > 0. 

(2) When g(A), the spectrum of A, is real, it is sufficient that xi > 0. 



(3) When a(A) c R, x, e (R p + )° © 0 q , and A = 





is reduced with 



p(An) > p(A 2 2), it is necessary and sufficient that A 22 e R q * q be quasimonotone 
nondecreasing with respect to a nonnegative cone. 

(4) When cr(A) cR,Xi e (r:>° © 0 q , and A is irreducible (or not reduced) it is 
necessary that for some permutation x, e (R p + )° © 0 q of Xi, and for the transformation 



B = 



( 1 ) 0 ‘ 



I 



with A) = B 'AS = 



X, a 
0 A 



, that the matrix A 22 be quasimonotone 



nondecreasing with respect to a nonnegative cone. In this case it is sufficient that X\ be a 
simple eigenvalue of A. 
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(5) When q(A) c R and A = 



satisfies the necessary conditions in (3) 



A n A 12 
0 A 22 

above, X = p(An) = p(A 2 2 ), and An and A 22 are in upper-triangular quasimonotone 
nondecreasing form with Ap+i ;P +i = X, then it is sufficient that A; iP n >0 for 1 < i < p. 

(6) For ct(A) cz R, it is necessary that Xj have positive components corresponding to Xi 
and at least all of the eigenvalues with nonnegative imaginary parts. Furthermore, 
considering just the irreducible quasimonotone forms (for which xi > 0) it is necessary 
and sufficient that the (shifted) spectrum of A solve the irreducible nonnegative inverse 
eigenvalue problem 

There are still some unresolved issues for the construction of the cone K(B). The 
case discussed in (5) above (Proposition IV.B.3) is not fully resolved. However, in order 
to apply the construction in (4) above (the proof of Theorem IV.B.2) to such matrices, the 
sufficient condition in (5) is also necessary. 

In Lemma IV. A. 1, constructing the matrix B _1 with nonnegative inverse B, using 
equation IV.A. la requires selecting a basis vector e^. Aside from trying all n 
possibilities, we do not have an efficient method for selecting e^. However, this is only 
needed when constructing the cone. The application of vector Lyapunov functionsdoes 
not require the cone’s construction, since it is sufficient to know such a cone exists to 
draw a stability conclusion. 

Furthermore, in Section IV. C we replace some e*, ke(l,...,p} with Xi when 
constructing the first cone. Again, we may have to try all p possibilities, although we did 
present a method of reducing the search. This is more important to applications that only 
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require the existence of a cone, as this constructive step is required before we can 
determine the quasimonotonicity of A22, and hence, A. 

A problem with the above technique, which is certainly a significant problem in 
applications that require constructing the cone, is that we must compute the generalized 
eigenvectors of a matrix. In case ( 4 ) above, if existence is all we have to determine, then 
it may be sufficient to compute the first two eigenvectors of A. However, to find the 
cone we need to compute all of the generalized left eigenvectors of A. While 
theoretically this construction always exists, and while our results give strong existence 
conclusions, it may be numerically intractible to actually construct the required cone. 

Concerning the related problems discussed in Chapter VI, for example the 
applications of M-matrices, in order to use the techniques of this paper to transform a 
given matrix into one of a particular form, then the cone K(B) must certainly be 
constructed. We do not yet know of a simpler construction than the one we presented. 

When the spectrum is complex we have the further problem of determining which 
part of the real spectrum is required to solve the irreducible nonnegative inverse 
eigenvalue problem Certainly we can use (for an irreducible, quasimonotone matrix) no 
more eigenvalues than there are nonzero components of xi, but we could have a further 
grouping of the complex eigenvalues with different collections of real eigenvalues when 
the final quasimonotone form of the matrix is reducible. We do not know a priori which 
collections may work, aside from trying all of them As in all of the above problems, since 
the method of vector Lyapunov functions is a proven method in examining large scale 
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dynamical systems (see Siljak, 1978; or Michel and Miller, 1977) we cannot expect to 
have simple numerical results when we employ these techniques. 

Finally, as long as the nonnegative inverse eigenvalue problem is unsolved, then 
this problem is unsolved for the complex spectrum as well, for we have shown that under 
the condition of a positive first eigenvector, the problems are equivalent. 

These questions present further directions for research, if not in terms of the actual 
solution to the problem, then at least in terms of the computability and efficiency of the 
method of solution, particularly when the application requires constructing the cone. We 
have provided analysts and researchers a basic theory and some useful techniques for 
approaching these problems. 
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